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The bulk properties and impurity effect in the one-dimensional S — 1/2 Heisenberg model with 
dimerization and next-nearest-neighbor exchange are studied by the density matrix renormaliza- 
tion group and exact diagonalization methods. This model widely covers the S = 1 Heisenberg 
model, two-leg ladders, bond-alternating chains and the 5* = 1/2 Heisenberg model. Almost 
the whole parameter space belongs to the Haldane phase. In the non-dimerizing limit a tran- 
sition to a kind of spin liquid takes place, characterized by the delocalized spin polarization. 
This transition causes the distinct impurity effects between a spin-Peierls compound CuGeOs 
and the Haldane systems. Experimental aspects of a bond-alternating compound CuNb2 06 are 
also discussed. Furthermore, the incommensurate correlation arising in this model is extensively 
studied, and a phase diagram is constructed. 
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§1. Introduction 

Recent progress of experimental studies has shed light 
on a variety of intriguing features in low-dimensional 
magnets. Above all, spin gaps found in such materials 
arouse broad interests with relation to the pseudo spin 
gap of the high-Tc superconductors. We list some classes 
of quasi one-dimensional (ID) spin-gap systems with typ- 
ical compounds: The Haldane systems (NENP), bond- 
alternating systems (CuNb206), spin-Peierls systems 
(CuGeOa), two-leg ladder systems (SrCu203), and so 
forth. Although these compounds have similar spin-gap 
properties without any magnetic order in low tempera- 
tures, the behaviors toward the elemental substitution 
are entirely different.^' ^^ A tiny amount of nonmagnetic 
or magnetic impurity (less than a few percent) breaks 
the spin gaps of CuGeOa^''') and SrCuaOg,^''^) and 
brings about antiferromagnetic (AF) long-range orders. 
In contrast, some Haldane systems"'^'^-' and CuNb206*-' 
keep their spin-gap behaviors even when substitutions 
as much as 20-40% are made, and no signs of magnetic 
orders are observed. 

As the origin of this difference, multiple factors can 
be considered: [1] Difference in the nature of spin states 
within chains, [2] higher-dimensional (interchain) mag- 
netic coupHngs,^'^"^ [3] couplings with lattice degrees of 
freedom, ^^^ and [4] other factors, for example, disorder 
induced by the impurities, etc. Among these, we notice 
the factor [1], allowing for the following facts. Various 
experimental results for the above compounds are sys- 
tematically and quantitatively understood within pure 
ID spin models.^' ^^-is) Furthermore, in spite of the sim- 
ilar interchain couplings (about 10% of the intrachain 
couplingi^'is)), CuGeOs and CuNb206 exhibit distinct 
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behaviors toward the impurity, as mentioned. 

The main purpose of this paper is to investigate the 
cause of the above impurity effect on the basis of ID 
chains. Using the density matrix renormalization group 
(DMRG), exact diagonalization (ED) methods and the 
perturbation theory, we study the properties within ID 
precisely, thereby discuss the interchain coupling and the 
lattice effect. In order to argue in a unified viewpoint, 
we consider a model which continuously connects all the 
above materials in its parameter space: The ID spin- 
1/2 Heisenberg model with next-nearest-neighbor (NNN) 
coupling and bond alternation, ^^"^^^ 
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(1.1) 



where J > 0, and a and /3 are dimensionless parameters 
controlling the coupling strength (see Fig. 1). 

First, we study the bulk properties of the ground state 
and low-lying excited states, to grasp the global feature 
the Hamiltonian possesses in its parameter space. Next, 
we study the properties of finite systems, namely seg- 
ments separated by the impurities, paying attention to 
the edge condition. In particular, the local spin polariza- 
tion is important to consider magnetic long-range orders. 
The key point of the present impurity issue is that in the 
transition area from the Haldane phase to the spin liquid 
phase the local spin polarization highly depends on the 
system size. 

Another fascinating aspect of this model is the incom- 
mensurate (IC) correlation. So far, except for the case 
without dimerization {j3 = 1)20-22) j^^ properties of the 
IC correlation have been elucidated. ^'^^ To complete the 
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phase diagram of the IC correlation in the a-[3 plane is 

another important purpose of this paper. 

This paper is organized as follows. In §2 we briefly 
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Fig. 1. Schematic diagram of the parameter space for the Haniil- 
tonian eq. (1.1), spanned by two dimensionless parameters: a is 
of the next-nearest-neighbor exchange and /3 of the bond alter- 
nation. 



review the previous studies on the model eq. (1.1). In 
§3 the bulk properties of several quantities are studied 
widely in the a-(3 plane. In §4 the incommensurate cor- 
relation is studied. In §5 we concentrate on the regime 
of spin-Peierls compounds (/3 ~ 1), and consider the im- 
purity effect in CuGeOa. In §6 the regime of Haldane 
phase is discussed with CuNb2 06 in mind. Finally, we 
summarize the main results in §7. In Appendix A, the 
relations between the two representations of the Hamilto- 
nian, eq. (1.1) and eq. (2.3), are given. In Appendix B we 
supplement the explanation in §3.1 on the anomalous be- 
haviors of the spin gap. In Appendix C we study the low- 
lying level structure for different boundary conditions, 
and show that the transition from the Haldane phase to 
the spin liquid occurs at /3 = 1. In Appendix D we study 
the dynamical properties of the bond-alternating chains, 
in addition to various properties of monoclinic CuNb2 06. 

§2. Model and Method 

The Hamiltonian eq. (1.1) is represented as a two-leg 
ladder with zigzag rungs as shown in Fig. 2. This model 
includes various classes of ID spin systems as described 
below. The parameter of bond alternation /3 is restricted 
to the range — cxo < /3 < 1, by taking account of the 
translational symmetry. For the open boundary condi- 
tion (OBC), however, it is quite important whether the 
first and last bonds have the value J or /3 J. This point 
will be discussed later. As for the parameter of the NNN 
exchange a, we restrict the range to < ct < oo in this 
paper. 

There are some special regions in the a-j3 plane as listed 
below (See also Fig. 1): 

[1] (a,/3) = (0, 1): The S = 1/2 AF Heisenberg model. 




Fig. 2. Connectivity of the Hamiltonian eq. (1.1). 



[2] /3 = -oo: The S' = 1 AF Heisenberg model. 
[3] /3 = 2a (0 < a < 0.5): The so-called Shastry- 
Sutherland (SS) line.^'*^ In particular, the two ends 
(0, 0) and (0.5, 1) are (a) the isolated dimer point, 
(b) the Majumdar- Ghosh (MG) point, ^^^ respec- 
tively. 
[4] /? = 0: The two-leg ladder. The isotropic ladder 

(Jieg/ Jrung = a — 1) is represented by (1,0). 
[5] a = oo: The model is decoupled to two S — 1/2 AF 

Heisenberg chains with the doubled unit cell. 
We start with the relation to the Haldane issue. On 
the (3 axis (a = 0), the Hamiltonian represents the AF- 
ferromagnetic (F) bond-alternating chains for /3 < 0. In 
the limit of /3 — > — oo (as the Hund coupling), TL is ex- 
actly reduced to the S = 1 Heisenberg model with a 
numeric factor 1/4. In this viewpoint Hida studied the 
properties for finite values of /? (a = 0), using ED.^^^ 
He found that the properties of the Haldane phase con- 
tinue even to the AF-AF alternating region (/? > 0). 
Later, by extending the non-local unitary transformation 
used for the S = 1 systems, ^^•' Kohmoto and Tasaki^^^ 
and Takada^^-' showed exactly that the hidden Z2 x Z2 
symmetry is completely broken for the case [3(a)], and 
approximately that the Haldane phase continues to the 
limit [1].3") 

On the other hand, there are many researches in view 
of the analogy between the Haldane systems and lad- 
ders. For example. White showed that the isotropic lad- 
der point [4] belongs to the same phase with the Hal- 
dane system, using DMRG.'^^-' Recently using a vari- 
ational matrix product method, Brehmer et al. stud- 
ied the ground-state properties^^^ and low-energy ex- 
citations^^-' of eq. (1.1), and concluded that their trial 
function smoothly connects the AKLT modeP^-' and the 
MG model [3(b)]. Thus, a wide range of the a-P space 
seems covered by the Haldane phase. In this phase the 
following conditions have to be satisfied: 

• A unique singlet ground state with a finite singlet- 
triplet excitation gap (Haldane gap) for the systems 
of L = 00 or for the periodic boundary condition 
(PBC). 

• Fourfold (pseudo-)degeneracy appears under the 
Haldane gap for the finite open chains due to the 
two localized S = 1/2 spins at both edges. ^'^-' 

• A hidden topological (string) long-range order 
(LRO) exists.34) 

At the MG point [3(b)], the ground state for the PBC 
spontaneously breaks the translational symmetry and 
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has twofold degeneracy: 



V'i-[1,2][3,4]...[L-1,L], 
V'2-[2,3][4,5]---[L,l], 



(2.1) 
(2.2) 



where [i,j] ~ a(i) f3(j) — (3{i)a{j) . On the other hand, on 
the SS Hne (3=2a (0 < a < 0.5), V'l whose singlet pairs 
sit on the bonds with the larger exchange J becomes the 
unique ground state, and -02 is no longer an eigenstate. 

On the a axis (/3 — l),^"-* the excitation is gapless for 
< a < ac. Note that the gapless regime in the a- (3 
plane is limited to this segment and a = oo. The criti- 
cal point Qfc = 0.241167^^) is connected to the Haldane- 
Shastry model, ^^•' in the sense that the conformal in- 
variance is realized. For a > ac, the ground state is 
twofold degenerate for L ^ oo, in accordance with the 
Lieb-Shultz-Mattis theorem. ^'^^ 

The IC correlation arises in the f3 < 2a side of the SS 
line.^°"^'^) However, details including the boundary posi- 
tion have been still obscure. As far as the a axis {(3=1) 
is concerned, it is known that an IC wave number which 
maximizes S{q)^^^ appears for 0.52063 < a < oo.^^'^^^ 

Here, we comment on another and more popular rep- 
resentation in the spin-Peierls regime {(3 ^ 1): 

i 

The Hamiltonians eqs. (1.1) and (2.3) are identical ex- 
cept that the space of eq. (2.3) is limited to 1 > /3 > — 1 in 
the a-f3 plane. We use each representation in its suitable 
situation. The relations between the two representations 
are summarized in Appendix A. 

Lastly, we describe numerical methods used. The 
DMRG method for zero temperature^^-* is a powerful 
method to calculate the ground-state and the low-energy 
properties accurately for low-dimensional systems. In 
this paper we adopt the finite size algorithm with OBC. 
The number of states kept for each block is up to 
m = 200. The system size is usually L ^ 200, and is 
increased up to 400 according to necessity. The typical 
truncation error is 0(10^^^) for the region away from 
the a axis (/? = 1). But it becomes 0(10^^) as /3 ap- 
proaches 1, where the correlation length diverges, and 
also in the incommensurate area. We show error bars 
when accuracy is required. 

The ED (Lanczos and Householder) methods arc used 
for supplementary calculations to DMRG, and to study 
dynamical and thermodynamic properties with a recur- 
sion procedure for the former. We check system-size de- 
pendence and make extrapolations from the data up to 
L = 26 for the former properties and up to L = 18 for 
the latter. 

§3. Global Features of the a-P Plane 

In this section, we study the bulk properties of several 
quantities widely in the a-/3 space. In §3.1 we consider 
the spin gap. The anomaly in the spin gap is a sign of 
the incommensurate correlation treated in §4. In §3.2 
two-point correlation function is discussed. In §3.3 the 
order parameters of dimcrization and of hidden topolog- 
ical (string) order are studied. We refer to the phase 



transition from the Haldanc phase to a spin liquid at 

/3=1. 

S.l Spin gap 

Before discussing the results, we need to remark on the 
edge effect in the DMRG calculations of the gap. As is 
known, an open 5* = 1 chain has free S = 1/2 spins near 
both chain ends in the Haldane phase. These spins yield 
pseudo fourfold degenerate states (one singlet and three 
triplets) beneath the bulk Haldane gap. Therefore, in 
estimating the bulk gap by DMRG,^"-* extra S = 1/2 spin 
had to be attached to both edges artificially to cancel 
the free spins. In the present case we take advantage 
of a similar but natural cancellation, by imposing the 
condition on the system that the edge bonds should be 
J (the stronger one) . For the bulk quantities, we use this 
boundary condition. We will discuss the effect of edge 
bond in Appendix C.2. 

Figure 3 shows the f3 dependence of the spin gap A 
for some values of a obtained by DMRG. The spin gap 
is obtained by A(L) = E{l;L)-E{0; L), where E{S''; L) 
is the lowest energy with total 5*^ and system size L. We 
extrapolate A(L) to the thermodynamic limit by using 
the polynomial fit: 



A(L) = A(^) + 5:|^, 



with n < 8 for i = 12, 16, 24, 32, 48, 64, 100, 200, 400. 

For a = 0, the result is consistent with those pre- 
viously estimated by ED.^^^ For /3 ^ — oo the gap is 
smoothly reduced to Ah/4J, where Ah = 0.41050(2)^^°) 
is so-called Haldane gap in unit of the exchange coupling 
for the S = 1 AF Hciscnberg chain. 

At the isolated dimcr point the gap becomes the maxi- 
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Fig. 3. /3-dependence of the spin gap for some values of a esti- 
mated by the DMRG method. The values are obtained by the 
extrapolation of the finite systems of L = 12-400. The inset 
shows the enlargement for 0.95 < /3 < 1- 



mum J in the a-(3 plane. There, the model is reduced to 
the two-spin problem, so that the energy cost to break 
up one singlet into triplet is J. At this point A forms 
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a cusp as seen in Fig. 3. In the following, wc study this 
anomaly in connection with the incommensurate correla- 
tion treated in §4. For finite but small values of a and /3, 
this anomaly can be traced by the second-order pertur- 
bation theory. The energy per bond for the ground state 
(singlet) and the lowest triplet branches are obtained re- 
spectively as, 



J 

EM 
J 



Es 3 3 2 

- = -4-T6(/3-2a). 



4 4^' 



Aa-" 



(3.1) 
(3.2) 



^•^ 



2a) + -(/32-4a2) 



cos(2<7). 



where q = mr/N {—N/2 < n < N/2; n = integer) and 
A'' = L/2 (the number of dimers) with the lattice con- 
stant being unity. Note that for (i = 2a cancellation oc- 
curs between a and [3. For the ground state, £^s leads to 
the constant exact value of the SS line.^'*^ Meanwhile, for 
Et{q) the hopping matrix element in eq. (3.2) vanishes, 
so that the singlet-triplet excitation forms a flat blanch. 
When a parameter [a or (3) is varied and crosses this 
line, the wave number minimizing Et{q) switches from 
q = tt/2 to (tt) discontinuously, due to the hopping 
term. Thus, the above dispersionless point corresponds 
to the cusp of the spin gap. On the other hand, the sec- 
ond term in eq. (3.2) survives even for /3 — 2a; the spin 
gap decreases from J even along the cusp line. 

When (a,/3) goes away from (0,0), the cusp still re- 
mains as seen in Fig. 3. In the a-(3 plane, the cusp point 
deviates to the side oi P < 2a with increasing a, as 
shown in Fig. 12. This singular point is related to the 
incommensurate correlation, and will be taken up again 
in §4. We have confirmed by ED that the properties 
found by the above perturbation theory are preserved 
for a ^ 0.5. On the other hand, for a ^ 0.5 with in- 
creasing (3, the bottom of the triplet blanches does not 
switch directly from tt/2 to (tt), but IC values of q inter- 
vene. Thus, the cusp disappears. Actually, for a = 1 the 
gap changes smoothly from the Haldane limit ((3 — — oo) 
to /3 ~ 1 through the isotropic ladder point {(3=0), as 
depicted in Fig. 3. The supplement of this issue is given 
in Appendix B. 

The inset of Fig. 3 shows the enlargement for 0.95 < 
P < 1. The (5-dependence^^' of the gap has been dis- 
cussed in relation to the spin-Peierls transition. ^■^''*^^ The 
gap drops steeply near (3 = 1 even for the spontaneously 
gapped regime {a > etc), as seen for a = 0.35 and 0.5. 

S.2 Spin-spin correlation 

In this subsection, we discuss the spin-spin correlation 
function Sij = {SfS^) in the ground state. Using the 
system of L = 280 typically, we calculate Sij in the cen- 
tral part of the system. The 40 sites near the open edges 
are not used to eliminate edge effects. 

In the gapped regime, where Sij decays exponentially 
for \i — i\ -^ oo,^^' we assume a simple form: 



{S^S^+^) oc exp - 



e 



Incidentally, we do not adopt the 2D Ornstein-Zernike 
type used for the S = 1 chain, ^'^^ because it is not certain 
whether it is adequate also in the present model. The 
above fit successfully works for most cases except for (3 ^ 
1, where the algebraic decay is essential. In Fig. 4 typical 
cases of the exponential decay are shown: An AF-AF 
and a F-AF bond-alternating systems and the isotropic 
ladder, etc. The inset shows that the ranges of ill fit 
by eq. (3.3) are limited to several lattice spacings. The 
correlation length S, is obtained from this fit. 
In Figs. 5(a) and 5(b) we plot £, for — 1 > /3 > — oo and 

10" 
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+ 
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(3.3) 



200 



Fig. 4. Spin-spin correlation function versus r for some typical 
cases of exponential decay. Since the fitted line for odd r shifts 
a little from that for even r for /3 > 0, only the data for even r 
are shown except for (0, —2.36). The inset is the magnification 
of the small-r regime with fitted lines. 



1 > /3 > — 1, respectively. In the Haldane limit {(3 -^ 
— oo), ^ is reduced to 2^h [cross symbol in Fig. 5(a)], 
where ^h — 6.03(1)^°^ is the correlation length of the 
5=1 AF Heisenberg chain. The factor 2 arises due to 
the doubled unit cell. As a increases, ^ approaches the 
value in the Haldane limit faster. This is because a does 
not work as frustration for /3 < 0, but plays a cooperative 
role to stabilize the ttii structure. 

As shown in Fig. 5(b) £, vanishes on the SS line 
{(3 = 2a), where Sij vanishes for |i — i| > 2.^^' Com- 
paring with Fig. 3, we find the minimum point of ^ does 
not coincide with the maximum point of the spin gap, 
except for the isolated dimer point. In the upper side 
of the SS fine (/3 < 2a), the region exists where \Sij\ 
weakly oscillates with an IC period. ^"^^ We indicate the 
other commensurate-IC boundary — so far unknown — by 
arrows in Fig. 5(b). In this regime, the fit by eq. (3.3) 
becomes less accurate with increasing (3 because of the 
frequent oscillation; possible errors thereof are shown by 
bars. We will discuss the IC correlation in §4. 

Now, let us consider the vicinity of the a axis. In 
Fig. 6 we show Sij for (3 = 1. For the gapless regime 
Sij shows algebraic decay (with some logarithmic cor- 
rection'*^)), as seen for a = and 0.2. An important fact 
is that even in the gapped regime {a > ac), Sij shows 
a similar non-exponential decay for some distance. Of 
course, in a sufficiently long distance Sij decays expo- 
nentially. Therefore, we try to divide Sij into two parts. 
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Fig. 5. Correlation length for some values of a as a function of 
(a) 1/0 for -co </?<-!, and of (b) /3 for -1 < /3 < 1. In (a) 
cross symbol at /3^^ = indicates twice the correlation length 
of the S = 1 AF Heisenberg model.'*"-' In (b) arrows stand for 
the commensurate-incommensurate boundary observed in Sij. ^ 
vanishes on the Shastry-Sutherland line (/3 = 2a). 



The arrows in Fig. 6 show the minimal r(= rj) which 
satisfies the condition that the error in the exponential 
fit for r > rj exceeds 10^'^ (this value is for convenience). 
Thereby, we roughly estimate the boundary of the decay- 
ing behavior. For a < ac, t] = oo; even for a > ac, t] is 
still considerably large, e.g. rj ~ 163 for a ~ 0.35. In the 
vicinity of the MG point {£, — 0) rj decreases abruptly. 

This tendency still remains in the vicinity of the a 
axis (/? ^ 1). There, the short-range correlation is con- 
trolled by the semi-algebraic decay, and poorly fitted by 
eq. (3.3). Therefore, the values of ^ plotted in Fig. 5(b) 
have been estimated from Sij^s with sufficient large |i— j|. 
Such enhancement of the short-range spin correlation 
near the a axis is closely connected to the impurity effect 
in spin-Peierls systems. This will be discussed in §5. 



S. 3 Order parameters 

We begin with an order parameter for dimcrization 
defined by 

Od = 1(52^-1 • ^2.) - {S2^ ■ S2^+l)\ ■ (3.4) 

Figure 7 shows Od versus [3 for some values of a. 




200 



Fig. 6. Spin-spin correlation function {SfSf.^) as a function of 
r for some values of a. /3 = 1. The arrow indicates rj, which 
roughly divides Sij into two parts, of exponential decay (r > rj) 
and of semi-algebraic decay (r < t]). The data for even r are 
shown for clarity. 



First, we look at the cases for a = 0. As (3 decreases 
from 1, the dimcrization abruptly becomes larger and 
continues increasing until (3 — —1, where Od has the 
maximum value 0.8245165. Since Ob = 3/4 for the iso- 
lated dimer point, moderate magnitude of ferromagnetic 
exchange between singlets further promotes the dimer- 
ization.^^) Beyond (3 = —\ Od monotonically decreases 
and tend to the value Od=0. 600371010 estimated for the 
S ~ 1 AF Heisenberg chain, for (3 -^ — oo. 

Next, for < a < ac there remain the absence of 
1.0 



Q 




Fig. 7. Dimer order parameter as a function of /3 (1 > /3 > —1) 
for some values of a. Inset shows the same quantity for /3 < 
— 1 as a function of inverse /3. The cross at /3~ = indicates 
the value for the S = 1 Heisenberg chain: (1 — Eii)/4, where 
E^ = -1.40148440390.40) On the SS line (/3 = 2a) On = 3/4. 
Extrapolation is carried out from the DMRG data of L = 12-400. 



dimcrization for (3—1 and the abrupt increase for /3 ~ 1. 
On the other hand, for a > ac^ Od shows a finite value 
even for (3 = 1 due to the spontaneous dimcrization. 
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With increasing a the niaxinium point of Od shifts in 
the positive direction of /3, although the maximum value 
is almost unchanged. This aspect is understood by the 
cooperative roles of a and — /3, as mentioned in §3.2. 
In contrast to the case of /3 = 1,^^-' there is no direct 
correspondence between the spin gap and Od. In Fig. 3 
the spin gap shows the cusp, while Od does not show any 
singular behavior, and smoothly connects the Haldane 
limit (/? = — cxo) and the limit of /? ^ 1. 

Now, we turn to the string order parameter Os- The 
string correlation function to measure the topological 
LRO in the Haldane phase for 5 = 1 systems^**) was 
extended to the S ~ 1/2 bond-alternating chain as,'^*'^ 



Os(*-.7) = 



2j-2 



(^2,-1 + ^2.) n ^M^^sf) iSI^_, + SIM. (3.5) 



e=2i+l 

In this expression, ferromagnetic (weaker AF) bonds are 
assumed between the sites 2i — 1 and 2i and between 
ones 2j — 1 and 2j. Among others, this definition is the 
most natural extension to the present system. 

In Fig. 8 shown is Os = linir-^oo Os('') estimated by 
DMRG. As has pointed out for the S* = 1 case,**"^ it is 
technically important to measure Os{r) at a sufficient 
distance from the open edge, not to mention large \r\. 
Actually, for the cases oi /3 ^ 1 and small a, the re- 
quired distance exceeds our chain length (L — 400). This 
technical difficulty prevents us from concluding decisively 
that Os vanishes as /3 ^ 1 for a < ac. We believe it oc- 
curs, however, taking account of the singular behaviors 
near j3 = 1. On the other hand, for a > a^^ Os has a 
finite value; for the MG point, Os = 1/4. 

When (3 decreases from 1, Os increases at first, has 
a maximum at a finite value of j3, and then tends to 
the value for the S* = 1 case. Here again, the value for 
a larger a approaches faster to the Haldane limit. For 
a — I, Os varies gradually through the isotropic ladder 
point (/3 = 0). Like Od, the string order changes con- 
tinuously keeping finite values from f3 = —oo to (3 ^ 1. 
This strongly suggests that the properties of the Haldane 
phase are kept in almost the whole a-(3 plane except just 
on the a axis. 

In fact, similar conclusions have been already reached 
in the case of a = by the ED and variational calcula- 
tions,^^' '^"' and for a ^ by variation methods. ^^'^^^ We 
demonstrate that the transition from the Haldane phase 
to a spin liquid phase occurs at /3 = 1 in the case of the 
SS line in Appendix. This discussion can be generalized 
to all the a-(3 plane. As we will see later, this "spin liquid 
phase" is characterized by the delocalized spin polariza- 
tion, in contrast to the localized moment at the edges in 
the Haldane phase. 



§4. Incommensurate Correlation 

In this section we study the incommensurate (IC) cor- 
relation the model eq. (1.1) presents. In §4.1 and §4.2 we 
discuss the IC correlation with respect to Sij and S{q), 
respectively. In §4.3 we construct a comprehensive phase 







Fig. 8. String order parameter as a function of /3 (1 > /3 > — 1) 
for some values of a. The inset shows the same quantity for 
— 1 > 13 > — oo. The cross at I3~^ = indicates O?(oo) = 
0.374325096(2) for the S = 1 Heisenberg model.""') On the SS 
line Os = 1/4. 



diagram of the IC correlation in the a-/3 plane. 

4.1 Disorder in the real space 

In classical spin systems with competing interactions 
which favor different wave numbers, e.g. n and 7r/2, a 
transition often occurs from an AF ordered state to spiral 
ones, as some parameter varies. Corresponding to this 
transition, in ID quantum spin systems the oscillatory 
period of Sij and the wave number which maximizes the 
static structure factor S{q) change from a commensurate 
(C) value to an IC one, although LRO's are not realized 
due to quantum fluctuations. Thus, incommensurability 
is measured in two ways, namely by Sij and S{q). Gen- 
erally, the transition point for Sij (disorder point) does 
not coincides with one for S{q) (Lifshitz point). ^ In this 
subsection, we focus on the former. 

Before discussing our results, we briefly review the 
classification and attribute of C-IC transitions. For 
Ising-spin systems, various properties of C-IC transi- 
tions were studied; the transitions are classified into two 
kinds. '*^' This classification is applicable to the quantum 
systems. ^^''*^-' The disorder point of the first kind has the 
following properties: In one side of this point \Sij\ has 
monotonic exponential decay, whereas in the other side 
the exponential decay is modified by an oscillatory factor 
with parameter-dependent wavelength, typically as. 



(5f5f+,)^ cos b(/3)r] £-'■/«(«. 



(4.1) 



In addition, the correlation length ^ behaves in the vicin- 
ity of this disorder point (/? = /3d) as. 



9C 



dp 



d^ 



dp 



finite, 



(4.2) 



IC 



where we assume /3 to be the parameter and C (IC) de- 
notes the limit of /3 ^ /3d in the C (IC) side. On the 
other hand, the disorder point of the second kind has the 
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properties that Sij decays exponentially with parameter- 
independent oscillations in both sides of the point. Be- 
sides, ^ satisfies the following conditions: 



e = (/3 = /3d), 



d(3 



(4.3) 



where the double sign indicates the limits from both sides 
of the disorder point. 

Now, we consider our model, which has the compet- 
ing AF couplings, the nearest-neighbor and NNN terms. 
As pointed out for Fig. 5(b), there are two C-IC bound- 
aries. The lower boundary coincides with the SS line. 
Hereafter, we use /3\ and /3JJ for the lower and upper 
boundary, respectively. As for the C areas, Sij has the 
ttJ-J- correlation (period 4) for /3 < /3J[, and the AF (stag- 
gered) correlation (period 2) for /3 > /3j. In Fig. 9, \Sij\ 
is shown for some typical IC cases (/3Jf < /? < /?j). In 
contrast to a smooth exponential decay for the C cases 
in Fig. 4, for /3 > /3JJ, \Sij\ has an oscillation with the 
long wave length which indicates the discrepancy from 
the period 4 (i.e. q = tt/2). With increasing /3, this os- 
cillatory period becomes shorter, as seen for a = 1 in 
Fig. 9. The inset shows the characteristic wave number 
q in eq. (4.1), obtained through the following fit: 

(5f 5f+,,)e«/'- ex cos [{q + n/2){r + 0.44)] . 

The exponent 9 oi q + tt/2 cx (/? — f3^y is estimated as 
9 = 0.39 for a = 1, for example. On the other hand, the 
condition for ^ eq.(4.2) is also satisfied, e.g. as seen obvi- 
ously for {a, (3) = (1,0.775) in Fig. 5(b). Thus, we have 
found the upper disorder line belongs to the first kind. 
Through the same procedure for (3\, we have confirmed 
that the lower boundary also belongs to the disorder 
point of the first kind. By the way, at the isolated dimer 
point the upper and lower boundaries merge into one, 
and the intermediate IC region vanishes [see Figs. 5(b) 
and 12]. Thus, in both sides of /3 = the oscillatory 
periods of Sij are constant (2 and 4), and eq. (4.3) is 
satisfied as seen in Fig. 5(b). Accordingly, the isolated 
dimer point belongs to the second kind. In this sense, 
the C-IC transition of the second kind is regarded as a 
complex of two transitions of the first kind. 

Finally, we summarize our results for three typical 
cases in Table I with those of the two related S = 1 
spin models for comparison. 



4.2 Incommensurate correlation in S{q) 

Another measure of incommensurability is the wave 



number q 
factor: 



which maximizes the static structure 



Siq) 



lE(^^^^l)' 



MiXi-Tj) 



(4.4) 



«j 



We call (7max as the pitch angle by analogy with the spiral 
spin order. The C-IC boundary (i.e. Lifshitz point) is 
defined as the point at which gmax deviates from the C 
wave numbers, tt or 7r/2. 

We have estimated S{q) by DMRG under OBC in the 



+ 
V 




Fig. 9. Spin-spin correlation function as a function of r for some 
parameter values in the IC area. Only the data of even r are plot- 
ted for clarity. The inset shows the characteristic q of eq. (4.1) 
for a = 1 (/3^ = 0.775). 



following way: Fixing the system size at Lq = 280, we 
calculate every Sij — {SfSp in the central part of the 
system. To eliminate the edge effect we do not use 40 
sites from both ends, so that S^'s of r^ = r^ — rj < 198 
are obtained. Additionally, thanks to the dimerization, 



Sij 's of negative r^ can be estimated only by shifting the 
starting site from i to i + 1. Therefore, we substantially 
obtain Sij for the system of L = 2 x 198. Since Sij decays 
exponentially in the region concerned, the long-distance 
contribution is negligible in the Fourier transformation. 
The values of S{q) thus obtained coincide with the ED 
results with PBC up to L = 26 within 7-8 digits. 

Figure 10 shows S{q) for various values of /?, and 
a = 0. In the 5 = 1/2 Heisenberg model (/3 = 1), S{q) 
is logarithmically divergent at q — n due to the quasi- 
long-range AF correlation.^"^ At the isolated dimer 
point (/3 = 0), S{q) = [1 - cos(g)]/4. For 1 > /? > 
/3l = —0.3723(2), (/max remains at tt, and starts to shift 
from TT at /3l (Lifshitz point). In the Haldane limit 
(/3 -^ —00), S{q) becomes a Lorentzian^^' sharply peak- 
ing at (/max = 7r/2. From the present calculations, how- 
ever, it is difficult to determine definitely the upper C-IC 
boundary. 

In Fig. 11 we plot (/max for some values of a versus (3. 
The pitch angle (/max suddenly drops from the higher C 
value TT, but slowly approaches the lower C value 7r/2, 
as /3 decreases. This aspect is common to the S = 1 
models.^*''*^) Since for a > 0.5206321) ^max is IC even 
for /3 = 1, (/max is always IC for o; = 1. In Fig. 11 also 
shown is (/max of the static susceptibility x{q) for a = 0, 
obtained by ED through, 

S{q,uj) 



x{q) 



duj- 



UJ 



This quantity is sometimes used to measure the incom- 
mensurability, but the pitch angle of x('Z) does not nec- 
essarily coincide with that of S{q). 

In this connection, the tendency of IC peak in S{q) 
has been observed in the inelastic neutron scattering 
spectrum for doped CuNb2 06,^^'' as discussed in Ap- 



Shinji Watanabe and Hisatoshi YOKOYAMA 



Table I. Comparison of the C-IC transition among the three models. The left column is the 5 = 1 bilinear-biquadratic Heisenberg 
chain: H = J [^ . S; ■ Si+i + 7^^(Si • Si+l)2] , where J > and 7 > O.-'^'**^) For 7 = 1/3 the model is so-called AKLT 
Hamiltonian^-^) and is written in the form of a projection operator. The middle column is the S = 1 frustrated Heisenberg chain: 
"H = J \ y^ , Si ■ Si+i + a y^ , Si ■ Si+i , where J > and a > 0.***' For the present model we tabulate three typical cases. In the 
section of Sij the values of the lower and upper disorder points are entered with their classes. In the section of S{q) the values of 
the lower and upper Lifshitz points are entered with their pitch angles [(jmax] in the C sides. In the 'Other transition' section we list 
transitions in the IC regime, if any. 



5=1 bilinear-biquadratic 
Heisenberg model 



5 = 1 frustrated 
Heisenberg model 



The present model 



a axis (/3 = 1) 



/3 axis (a = 0) 



a = 0.35 



upper 



7 = 1/3 (AKLT) 

(first kind) 

7 = 1 



a = 0.284(1) 
(first kind) 



5^ = 0.5 (MG) 
(first kind) 



/3 = 

(second kind) 

13 = 



/3 = 0.7 

(first kind) 

fS = 0.445(5) 

(first kind) 



5(g) lower 
upper 


■y = 0.43806(4) [vr] 

7 = 1 [27r/3] 


a = 0.373(3) [n] 
a = 00? [7r/2] 


a = 0.52063(6) [n] 
a = 00 [7r/2] 


13 


= -0.3723(2) [tt] 
/3 = 00 ? [7r/2] 


13 = 0.5131(1) [tt] 
/3 = 00 ? [7r/2] 


A(max) 
Os(max) 


7 = 0.406 
7 = 1/3 


a = 0.40(1) 
a = 0.4397(1) 


a~ 0.6 
a = 0.748(2) 




/3 = 
/9 = -4.47(1) 


13 = 0.525(5) 
/3 = -2.93(1) 
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Transition 


No 
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Fig. 10. Static structure factor for various values of /3 and a = 0. 
The system used corresponds to the periodic chain of L = 396 
(see text). Note that the line only connects the discrete data 
points without any interpolation. 



Fig. 11. Pitch angle gmax/Tr of S{q) as a function of /3 for several 
values of a. Since the data are of L = 396, only discrete values 
of Qmax are allowed. Filled circles shows gmax/Tr of the static 
susceptibility x{l) for a = 0. 



pendix D. 

4-3 Phase diagram of incommensurate correlation 

So far, we have studied the IC correlations in terms of 
three quantities: The wave number of the lowest triplet 
excitation Et{q) in §3.1, the spin-spin correlation Sij in 
§4.1 and the static structure factor S{q) in §4.2. For each 
quantity we have determined C-IC boundaries in the a-f3 
plane, which are depicted in Fig. 12. 

The dashed line with triangle represents the boundary 
between q = tt (0) and tt/2 for the lowest Et{q). On this 
line the spin gap shows the cusp as seen in Fig. 3. This 
jump from tt (0) to 7r/2 with decreasing /?, i.e. cusp, is 
obvious for a ^ 0.5, but it disappears for a ^ 0.5 due to 
the IC excitations (see Fig. 24 in Appendix). This line is 
included in the IC region of Sij , and discords from each 



boundary mentioned below. 

The shaded area is the IC region where | Sij \ decays ex- 
ponentially with a quasi-periodic oscillation as in Fig 9. 
This region is bounded by the disorder lines indicated 
by the solid lines with filled circle. The C region un- 
der the lower (over the upper) disorder line has the AF 
(ttii) correlation. On the lower boundary (the SS line), 
the Hamiltonian can be written in the form of projection 
operators. '^^•' The fact that Ti, is represented by projec- 
tion operators at the disorder point is common to some 
other models. In the spin-orbit model the disorder point 
is the SU(4) point, ^^^ and in the S" = 1 model the AKLT 
point, *®^ as shown in Table I. 

The thin solid lines show the contours of the pitch an- 
gle 9max/7r, which maximizes S{q). In contrast to other 
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cases, '*^''**^ the IC region of S{q) extends widely in the 
C region of Sij . The pitch angle decreases abruptly near 
the Lifshitz line (gmax = tt). In such a region S{q) shows 
merely a broad hump instead of a sharp peak (Fig. 10), 
hence ^max is hardly regarded as a characteristic wave 
number there. Indeed, we have not observed a singular 
behavior in any quantity at the Lifshitz point; system- 
size dependence of many quantities in this IC region 
shows the period of 4, in accordance with the C period 
of Sij. Thus, the intrinsic IC correlation is represented 
in terms of Sij . 

In Fig. 12 all the lines go down rightward. As repeat- 
edly mentioned, this indicates that the properties are 
preserved along this direction due to the similar role of 
—f3 and a. 

Lastly, we touch on the continuity of the Haldane 
phase. The IC regime lies between the Haldane limit 
(P = — oo) and /3 = 1 irrespective of the value of a. In- 
deed, the correlation length and the spin gap are singular 
there. However, various quantities including the dimer 
and string order parameters (§3.3), which characterize 
the Haldane phase, have no anomaly. Actually, we will 
confirm that the Haldane phase extends to the limit of 
/? ^- 1 in §5 and Appendix C. Hence, the IC correlation 
does not affect the essence of the Haldane phase and also 
of the spin liquid phase. 



cusp of spin gap A 
incommensurate boundary 
determined by <SfS/> 



ImJ^ofSiq) 




Fig. 12. Comprehensive phase diagram of IC correlation in 
eq. (1.1). Shaded area is the IC regime in Sij. Thin solid lines 
are the contours of pitch angle gmax which maximizes S{q); the 
Lifshitz line corresponds to the line of gmax/7r = 1.0. The IC 
area in S{q) extends further in the directions of o — > oo and 
P -^ — oo. Dashed line with open triangle shows the cusp of the 
spin gap, which distinguishes the wave numbers of the lowest- 
energy excitation between q = n and 7r/2 (see §3.1). 



§5. Spin-Peierls Regime 

In this section, we concentrate on the Spin-Peierls 
regime. First, the low-lying excitation is reconsidered 
in §5.1. Next, we study the edge effects in the light of 
the delocalized spin polarization, and consider the origin 
of the impurity effect of CuGeOs in §5.2. The stabil- 
ity of the edges is considered in reference to the lattice 
distortion in spin-Peierls compotmds in §5.3. The NMR 
spectra is discussed in §5.4. 

5.1 Two-spinon bound .state 

One primitive modeling of the substitution by non- 
magnetic elements is to consider finite-size systems, as- 
suming that the impurities simply cut the ID chain into 
finite segments. For finite systems, the boundary con- 
ditions are essentially important, in particular for the 
dimerized cases (/3 < 1). Then, there are two kinds of 
exchange coupling at the chain ends: J (stronger one) or 
PJ (weaker one). For simplicity, we call the former (lat- 
ter) case as the strong (weak) edge. In even-site systems 
both of the edge couplings are J or (3 J. 

The local spin polarization (S^) is an important quan- 
tity to study the local properties of magnetic excitation 
and the tendency toward magnetic LRO's.^^-' Figure 13 
shows (5f ) of the lowest-energy state in the S^ = —1 
subspace for a = (I) and a — 0.35 (II). For both values 
of a the cases of the strong and weak edge arc shown in 
(a) and (b), respectively. As described in Appendix C.2, 
the strong edge has some analogous properties to PBC 
or to the bulk. Since the main subject of this subsection 
is the low-lying excitation in the bulk, we restrict our 
discussions to the qualitative feature of the strong edge, 
here. We will concentrate on the edge properties and 
impurity effect in §5. 2- §5. 4. 

In Figs. 13 1(a) and 11(a) there are two humps of stag- 
gered moment, and the amplitude becomes maximum 
for P = 1. These humps are considered as the amplitude 
of excited free spinous {S — 1/2 entity). ^"^'^^^ For the 
S = 1/2 Heisenberg model [(a,/3) =(0,1)], it is known 
that the excitation can be exactly described by the scat- 
tering of even number spinons.^^' ^*-' For the MG point 
[(a,/3)=(0.5,l)], a similar free spinon picture was suc- 
cessfully introduced from a variational viewpoint,^'*-' and 
recently it was numerically shown that the MG gap (two- 
spinon excitation) is accurately twice the bottom of the 
single spinon blanch for odd-site systems. ^■^•' Thus, the 
free-spinon picture seems applicable to all the values of 
a for /3 = 1 at least in the low-lying excitation. When 
the dimerization is switched on (/? < 1), the amplitude 
of the oscillating part is abruptly reduced, and almost 
disappears for P = 0.2 leaving a broad peak as seen 
in Figs. 13 1(a) and 11(a). The mechanism of this pro- 
cess has been pursued by many studies.^*' ^^'^^"®^'' When 
P decreases from 1, two spinons are bound to make a 
magnon (5 = 1 entity) described by a single parameter; 
an isolated blanch thereby appears below the continuum. 

Furthermore, multi-magnon branches have been found 
above it. "'53, 62) 

To grasp intuitively the process from spinon to 
magnon, it is useful to consider the excitation along the 
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Fig. 13. Local spin polarization of the lowest-energy state in the subspace of 5"^ = — 1 for open systems with L=200 as a function of 
i for a=0 (I) and q=0.35 (II). At both chain ends the strong coupling J is assumed in (a), while the weak coupling f^J in (b). For 
/3 ~ 1 in every panel, {S^) oscillates alternately with considerable amplitude. 



SS line. At the MG point the ground state is twofold 
degenerate, ipi and ■02 in eq. (2.1). Here, we call the 
singlet array of V'l (V'2) as the R (L) array [Figs. 14 1(a) 
and 1(b)]. In the lowest-energy excitation a singlet pair 
breaks into two spins (spinons), which sit on the mutu- 
ally different subchains (legs) to avoid breaking another 
singlet pair, as shown in Fig. 14 II. The resultant spinons 
become two domain walls, between which the phase of 
singlet array, say L, is opposite to the original one, R. In 
this case the L array has an equivalent energy to the R 
array, so that the spinons move freely regardless of the 
distance between them. Thus, these scattering states 
make a low-energy continuum. On the other hand, when 
P is reduced from 1 keeping /3 = 2a, the degeneracy of 
the ground state is lifted. One of the states, say V'l in 
which singlets are positioned on the strong bonds (J) is 
still the exact ground state for /3 < 1, while ip2 is no 
longer an eigenstate and even has high energy over the 
singlet-triplet gap, as explained in Appendix C. Con- 
sequently, the L array is energetically unfavorable and 
works as a linear confinement potential; the two spinons 
approach each other to reduce the region of L. Namely, 
the configuration of Fig. 14 Ill(b) is more likely than that 
of Fig. 14 Ill(a). The resultant bound state causes an 
isolated branch below the continuum. Although we have 
focused on the SS line, this explanation is qualitatively 
valid for all the other paths from the a axis. This change 
of the excitation mechanism corresponds to the transi- 



tion from the spin liquid phase to the Haldane phase, ^*' 
and is analogous to the so-called A chain. ®^' 

Actually, one can confirm the appearance of the isolated 
magnon branch in the dynamical structure factor, 

S{q,i0) = Y, \{'^n\S'g\'^o)f5 {lo - {En - Eo)) , (5.1) 
n 

where h = I, ^n (En) is the n-th eigenfunction (eigen- 
value) of the system, and the ground state is indicated 
by n = 0. n runs over all the eigenstates. S{q,uj) is 
proportional to the intensity of inelastic neutron scat- 
tering spectra. In Figs. 15(a)-(c) we show S{q,uj) for 
a — 0.35, which is the plausible value for CuGcOs.^'*-' As 
discussed in the previous paper, ^'^•' the lowest mode for 

5 = Q {fi = I) [Fig. 15(a)] is assigned to the contimuim 
by checking the system-size dependence of the residue 
and position of pole in S{q,io). As soon as the dimer- 
ization is introduced, weak intensity appears below the 
continuum, whose position is almost independent of L, as 
seen in Fig. 15(b). This branch remarkably grows when 

6 increases to 0.012 (/3 = 0.98) [Fig. 15(c)], which is an 
expected value for CuGeOs within eq. (1.1).^^'^^^ This 
abrupt increase of the intensity is obvious in the inset 
of Fig. 15(b), where the ratio of intensity (residue) for 
the lowest branch at q = 7r/2 is plotted. Since the inten- 
sity belonging to the continuum decreases with 1/L as 
L ^ 00, this isolated magnon branch makes a dominant 
contribution. 
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Fig. 15. Dynamical structure factor S{q,uj) for three values of <5 are shown for the fixed value of o = 0.35. Data obtained by the ED 
with a recursion method for L = 14-26 are simultaneously plotted. The intensity is proportional to the area of circle. The magnon 
branch is indicated by arrow in (b) and (c) . The inset in (b) shows the intensity ratio (residue) of the lowest branch at q = it/2 versus 
5 for L = 20 and 24. System-size dependence is small in this wave number. 



Fig. 14. Schematic figures for low-energy excitations along the 
SS line (/3 = 2a). 1(a), (b): The degenerate ground states for 
the MG point. II: The lowest-excited state at the MG point, 
in which two spinous on different subchains (legs) form domain 
walls. III(a),(b): Excited states on the SS line (/3 < 1), where 
two spinous are located apart and closely, respectively. See text 
for details. 



Recent neutron expcrinients^®-* for CuGe03 vividly 
showed the appearance of the magnon branch. In the 
measurement at T = 12K (immediately below Tsp ^ 



14K) the spectrum has a sharp peak below the contin- 
uum as previously observed,^''' whereas at T = 16K this 
sharp peak suddenly vanishes and only a broad contin- 
uum is observed. Thus, the change in the excitation 
mechanism between 6 — and S > has been confirmed 
by experiment. 

Lastly, we touch on the thermodynamics. According 
to the abrupt change of the level structure near 6 = 0, 
thermodynamic quantities drastically changes, too. We 
confirmed that by calculating the specific heat C{T) 
along the SS line. As soon as a decreases from 0.5, a 
sharp peak appears. This peak is attributed to the dense 
low-lying triplet levels caused by the dimerization. The 
situation is the same with the A chain. ^'^•' 

5.2 Delocalized spin polarization 

In this subsection, we study (S^) in the lowest mag- 
netic state (5 7^ 0) to clarify the impurity effect. We as- 
sume that non-magnetic impurities cut the chain at ran- 
dom into odd- and even-site segments. This assumption 
is not so crude, since the NNN exchange hardly change 
the situation, as will be discussed later. 

First, we look at the odd-site systems. Although for 
even-site segments the ground state is always singlet and 
(Sf) = for each i, yet for odd-site segments the ground 
state is S* = 1/2, therefore (5f) ^ 0. In Fig. 16(a) we 
show {S^) for L = 21 and a = 0.35 in the S'- = 1/2 
sector. For the non-dimerizing case (/3 = 1), the spin 
polarization extends to the whole system, corresponding 
to a spinon discussed in §5.1. This delocalized spin po- 
larization is still seen in the IC regime (a > 0.5) as far 
as /? — 1, as shown in Fig. 16(b). Therefore the delocal- 
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ized spin polarization is characteristic of the spin-liquid 
phase, and has nothing to do with the existence of a 
gap. For /? < 1 one of the chain ends is the strong edge 
[i — 21-22 in Fig. 16(a)] and the other the weak edge 
(i = 1-2). With decreasing (3, the magnetic moment be- 
comes localized near the weak edge and the amplitude 
near the strong edge vanishes. Here, let us return to 
Figs. 13 1(b) and 11(b), where the lowest excited states 
(S — 1) in even-site segments with the weak edges are 
shown. These states are obviously regarded as the dupli- 
cation of the ground state of odd-site systems (S = 1/2) 
connected at the strong edges. Hence, the localization 
of magnetic polarization is a characteristic of the weak 
edge. This free 5* = 1/2 excitation at the edge is nothing 
but the celebrated aspect of the Haldane chain. ^'^'^*) 

In the following, we treat this localized moment more 
quantitatively. Consider a finite system with L spins, as 
in Figs. 13 1(b) and H(b). For the non-dimerized case 
{P = 1), the range of finite {S^) extends to the whole 
chain irrespective of the value of a, as mentioned. When 
dimerization is introduced, this range shrinks toward the 
edges. Incidentally, in the spin-Peierls regime (/3 ^ 1), 
(Sf) around the edges does not obey the simple expo- 
nential decay unlike the Haldane regime, so that the lo- 
calization length is not an appropriate measure of the 
localization range. Therefore, we introduce a somewhat 
artificial quantity r/odgc, defined as the maximum inte- 
ger satisfying KS"?)] > 10^"^ for all i < ryodgo- Here, we 
choose the value 10""^ merely for convenience. In Fig. 17 
we show 77odge estimated for the systems of L = 200. 
Note that ?7odgc scarcely depends on the system size, in 
contrast to (S'f ) for the strong edge. The localization 
range abruptly becomes narrow with decreasing (3 from 
1. With increasing a, ?7cdgo behaves more abruptly. For 
/3 ^ 0.9, (Sf) around the edges shows the typical expo- 
nential decay characteristic of the Haldane phase. 

Next, returning to the Figs. 13 1(a) and 11(a), we give 
quantitative treatment of the strong edge. In this case 
the local moment can be decomposed into the uniform 
and the staggered parts as,^^' '''^' 



{s!)^{s!r + i-iy{s!r 



(5.2) 



For a = the uniform part (S^)^ is almost constant es- 
pecially in the bulk for /3 = 1; (S'f )" abruptly deforms 
with decreasing /?, and converges to a sinusoidal curve for 
/3 -^ 0. For instance, see /3 — 0.2 in Figs. 13 1(a). Since 
the amplitude of the curve is determined by the magne- 
tization and the system size, (S'f )" oc 1/L as shown in 
the inset of Fig. 17. On the other hand, the staggered 
part {S^y is drastically suppressed with decreasing /3 
and vanishes at /3 = 0, as shown in Fig. 17. Further- 
more, it is important that (Sf)^ highly depends on the 
system size, and vanishes as L — > oo (inset of Fig. 17). 
Thus, {Sf) - (S'f )" for /3 - 0. The behaviors of (Sf) is 
basically the same for a = 0.35. 

Based on the above results, we consider the impu- 
rity effects in the spin-Peierls regime, especially with 
CuGeOs [a = 0.35, /? = 0.97-0.98 for T < Tgp] in mind. 
For pure CuGeOa a magnetic LRO has never been ob- 
served, while for a doped sample with 0.12% of Zn shows 



Fig. 16. Local spin polarization in the ground state (5^ = 1/2) 
for various cases, (a) Segment with odd sites for some values 
of fd. A weak edge (/3J) is assumed in the left chain end. (b) 
Segment with L = 21 in the IC regime without the coupling 
alternation, (c) System with L = 24 in which the 12th site is 
replaced by a site impurity. See text for details. In the margin of 
(a) and (c), the respective connectivity is schematically shown. 
As for (b), the connectivity is the same with (a). 



Fig. 17. Amplitude (maximum value) of (S^)^ for the strong 
edge, and the range of localized moment near the chain ends 
^cdgo for the weak edge, as a function of /3. The inset shows the 
system-size dependence of the uniform part (value at the mid- 
chain site) and the staggered part (maximum value) of (S^) for 
the Hcisenberg model. The used states have the lowest energies 
in the 5"^ = — 1 subspace. 



AFLRO in the low temperature.^-' 

[1] Infinite systems (L = oo, pure material) belong to 
the Haldane phase, as discussed in Appendix C. (a) For 
the open strong edge, the staggered polarization extends 
to the whole chain, however its amplitude becomes in- 
finitesimally small as L ^ oo, as shown in Fig. 17. Thus 
a magnetic LRO will hardly arise, even if the interchain 
coupling is somewhat large, (b) For the open weak edge, 
there appears staggered polarization with considerably 
large amplitude in the range of about 30 sites around 
the edges, but this range is negligible against L ~ oo. 
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The bulk area is not affected by this moment; the higher 
dimensional coupling will hardly induce a magnetic LRO 
also in this case. Thus, regardless of the kind of edge, 
infinite (pure) systems will keep the gap properties of the 
Haldane phase (see §6). 

[2] Finite systems (impurity-doped material) may 
break the characteristics of the Haldane phase (Appendix 
C). (a) For the open strong edge, the amplitude of (5f )'' 
in the whole system increases with the system size de- 
creasing, in particular abruptly for the thin impurity con- 
centration (inset of Fig. 17). Thus, making use of the en- 
hanced staggered moment, the interchain coupling prob- 
ably induce AFLRO. (b) For the open weak edge, with 
decreasing the system size, the bulk part with (S'f) = 
vanishes, since ryodge is independent of L. As a result, 
the localized polarization around the edges extends to 
the whole system. AFLRO is likely to appear also in 
this case. 

To summarize, when the system becomes finite by 
doped with non-magnetic impurities, the properties of 
the spin liquid are brought about to the spin-Peierls 
regime, which originally belongs to the Haldane phase 
for the pure case, and AFLRO will be induced by the 
interchain magnetic coupling. The above results are not 
qualitatively influenced by the value of a, so that the 
same mechanism works for spin-Peierls materials in gen- 
eral. So far, we have not made allowance for the stabil- 
ity of the edges. In fact, the strong edge is always stable 
against the weak edge in the spin-Peierls regime. We will 
discuss this point in the next subsection. 

In this connection, recent /iSR measurements for 
doped CuGeOa have observed the spatial inhomogeneity 
of the ordered moment size.^®-* However, its local struc- 
ture has not been elucidated yet. Meanwhile, concerning 
(Sf), the result of the phase Hamiltonian approach™^ is 
similar to the weak-edge case, if anything. 

Finally, we mention the effect of the NNN coupling 
across the impurity. Consider a system consisting of 
24 sites but the spin on the 12th site is replaced by a 
site impurity like Zn in CuGeOa [Fig. 16(c)]. Thereby, 
the system is separated into odd-site (11) and even-site 
(12) segments by the impurity, but the NNN exchange 
{a = 0.35) connects them. Since this system has odd 
spins (23), the ground state is magnetic {S = 1/2). In 
Fig. 16(c) we shows a typical result among various cou- 
pling conditions calculated. Here, we set the weak edge 
in the odd-site segment next to the impurity (between 
i = 10 and 11) and the strong edges for the even-site 
segment. An important aspect of this result is that the 
large amplitude of (S^) in the odd-site segment scarcely 
extends to the even-site segment beyond the impurity. 
This aspect does not change under other coupling con- 
ditions and also in the cases of a bond impurity like Si 
in CuGeOs. Summing up, the NNN coupling does not 
severely affect our simple treatment of the impurity ef- 
fect. 

5. 3 Stability of the edges 

An important factor for the impurity effect in spin- 
Peierls compounds is the softness of the lattice. In the 
low temperature the lattice distorts so as to minimize 



the sum of the lattice distortion energy and magnetic 
one. We can allow for this factor to some extent by 
considering the stability among different edge conditions. 
First, we consider the energy difference between the 
ground states of the strong and the weak edges in even- 
site systems. In Fig. 18(a) we plot 



AS(L) = 



.B(SS,L)-S(WW,L), cvenL 
.B(SS,i)-i;(SW,L), oddL 



(5.3) 



versus 5 for a = 0.35. Here, i<^(SS, L) means the ground- 
state energy of the system with L sites and the strong 
edges [cf. Fig. 18(b)]. In this range of 5 the energy of 
the strong edges is always lower than that of the weak 
edges. The behaviors of £;(SS, L) and £;(WW, L) are 
identical with those along the SS line. Namely, i?(SS, L) 
is a monotonically decreasing function of 5 and propor- 
tional to L, whereas _E(WW, L) highly depends on L 
for (5^0, and has discontinuity a,t 5 — Q for L ^ oo. 
Therefore, AE{oo) in the limit of (5 ^ is finite and 
equivalent to the gap at (5 = 0. For details, refer to 
Appendix C. The finite A_E(oo) for (5 — > in this case 
is recognized in the inset, where dAE/dS diverges for 
L ^ oo. Next, for odd L we also plot AE{L) of eq. (5.3) 
in Fig. 18(a). Here, i?(SW, L) is the ground-state energy 
of regular odd chain, while _E(SS, L) is that of the strong 
edges at both chain ends [See Fig. 18(b)]. The defect (ir- 
regular pattern) of exchange coupling arising in the latter 
is placed at (or near) the center of the chain. Since AE 
is negative even in this case, the strong edge lowers the 
energy, overcoming the energy loss by the defect. Inci- 
dentally, -E(WW, L) is much higher than both energies 
above. Furthermore, concerning the first excited (triplet) 
state, we compare the various cases of Figs. 13(a) and (b) 
{L ~ 200). Also in these cases the strong edge has lower 
energy than the weak edge. 

Thus, we conclude that in the spin-Peierls regime both 
energies of the ground state and of the lowest excited 
state are stabilized by the strong edges. In other words, 
if the lattice is sufficiently flexible like spin-Peierls com- 
pounds, the edge becomes strong. We would like to em- 
phasize that this conclusion holds even when the elas- 
tic energy of the 3D lattice is included. Because, as 
is known, the decrement of magnetic energy due to the 
spin-Peierls distortion E'mag oc (5" (j^ < 4/3) overcomes 
the increment of elastic energy of the lattice i?iat oc 5"^ 
for ,5 -^0.13.71) 

The result that the strong edge is favorable is consis- 
tent with the enhancement of singlet correlation at the 
edge bond,"'''''^'^) although this is not the case for large 
negative values of /3, as will be considered in §6. In the 
above we have discussed simply with flxed-frame lattices. 
To treat the spin-Peierls transition more accurately, one 
needs to introduce the magneto-elastic coupling. How- 
ever, the results of recent studies^^' allowing for such a 
factor are basically the same with ours. 

Based on the above results, we speculate on the first- 
order compositional phase transition in CuGeOa in the 
remainder of this subsection. Recently, Masuda et alJ^^ 
reported that when non-magnetic impurity Mg is doped 
as Cui_a;Mga;Ge03, a clear first-order transition is ob- 
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Fig. 18. (a) Energy difference of the ground state between the 
two edge conditions, which are different between the even and 
odd segments, as shown schematically in (b). For details see text. 
In the inset of (a) plotted is the inverse of JdAE/dS estimated 
by small values of 5. 



served at x = Xc ^ 0.023 (later modified as 0.027^^^) 
between the dimerized and uniform AF ordered pliases. 
Tlie dependence of Tn on x has a discontinuity at Xc and 
the spin-Peierls order disappears for x > Xc- A recent 
neutron scattering study^"'' has supported this first-order 
transition. According to it, the amphtude of magnetic 
moment fx^s as well as that of the spin-Peierls distortion 
have jumps at Xc- 

For T < Tsp the lattice of CuGeOs distorts in every 
direction of the three crystal axes^**) (The c axis is the 
chain direction) . The array of bond alternation in a chain 
is "out of phase" against those in the adjacent chains in 
both a- and b-axis directions, namely -L-R-L-R- in the 
notation of Fig. 14 I. When a non-magnetic impurity 
cuts the chain, the two edge bonds next to the impurity 
have to be strong. Thus, the array in one side of the 
impurity becomes "in phase" locally, e.g. -L-L-L-. This 
"in-phase" array raises the elastic energy, so that this 
reversed phase should be relaxed to the original "out- 
of-phase" array in some distance from the impurity.^^' 
Actually, it is observed that the width of spin-Peierls su- 
perlattice peaks in the neutron scattering spectra shows 
the resolution limit, ^^•' namely the intrachain phase of 
dimerization docs not change globally across the impu- 
rity. 

When the concentration of impurity is thin, the con- 
tribution of the above local defects to the total energy is 
smaller than the energy gain by the global spin-Peierls 
distortion. As x increases, local defects appear here and 
there, which damage the latter energy gain. Eventually, 
for X > Xc the contribution of the defects overcomes 
the energy reduction by the global dimerization. Thus, 
the compositional transition of the first order may be 



understood by the competition between the 3D global 
spin-Peierls distortion and the ID local defects. ^^•' The 
discontinuous increase of HeS at a; = Xc is understood 
naturally; if the dimerization vanishes by the above first- 
order transition, (S^) is enhanced discontinuously. For 
instance, imagine that /3 jumps from 0.99 to 1.0 in Fig. 13 
11(a). 

5.4 NMR spectra 

In this subsection we discuss the impurity effect on the 
NMR spectrum g{H). Since non-magnetic states [S = 0) 
do not contribute to the local susceptibility, 

1 



Xi 



T ^-^ 

3 



{SiS^)T, 



(5.4) 



Xi is proportional to (5f ) with respect to the lowest mag- 
netic state in the T ^ limit. In eq. (5.4) (• • ■)t indi- 
cates the thermal average. Then, by replacing Xi{T) by 
(Sf) the NMR spectra is expressed as 



g{H) = Y,f[H-{St) 



(5.5) 



where f{H) is the NMR line shape with natural line 
width at each site, and we assume f{H) = exp{—H'^/w^) 
in this paper. And H — H/{Ag^B) with A being the 
hyperfine coupling of the material. Like eq. (5.2), local 
susceptibility is sometimes decomposed into the uniform 
and staggered parts: 



x.-x-+(-irx- 



(5.6) 



When xf is independent of z, it is nothing but the bulk 
susceptibility x — xf ^i^d only shifts g{H). Actually, 
X" often depends weakly on i as in Fig. 13(a). On the 
other hand, x\ determines the line shape. When x\ has a 
different oscillatory period from tt, e.g. for the IC or ttii 
structure, the decomposition of eq. (5.6) is unavailable. 
First, we look at the Heisenberg model (a = 0, /? = 1). 
Eggert and Affleck have obtained Xi{T) for OBC by the 
conformal field theory i^^-* 

X%T) = -^ (5.7) 

TT y(j72Tyih5ri4TlA7) 

with a constant a ^ 0.58. The chain end is at i = 0. 
In relatively high temperatures Xi has finite amplitude 
only near the chain end, whereas with decreasing tem- 
perature the oscillation extends to the inner part of the 
chain and its amplitude increases. This aspect is con- 
sistent with the delocalization of {Sf) for T ^ and 
(3—1. Those changes cause the broadenings of NMR 
line shape and of background intensity with sharp edges, 
respectively. Such behaviors were actually observed by 
an NMR experiment for a good Heisenberg compound 
SrzCuOa.^s) 

Since eq. (5.7) is the result for the infinite system, we 
consider finite-size effects. In Figs. 19 I and II we show 
g{H) for a — and 0.35 respectively, estimated by the 
data in Fig. 13 I and II. For /3 ~ 1 and regardless of 
the value of a, g{H) has the double-peak structure on 
account of the delocalization of {Sf) with mild curva- 
ture. Indeed, such a structure has been observed in the 
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NMR spectrum for Sr2Cu03 at a low temperature [see 
Fig. 2 (30K) of rcf. 78]. In this spectrum the interval 
between the peaks looks narrower than that in Fig. 19 I 
(L = 200). This is interpreted as follows. The amplitude 
of {S^y (and also (•S'f)") becomes smaller and flatter 
with increasing L, as explained in §5.2. Thus, the inter- 
val becomes narrower with increasing L, as observed in 
Fig. 20(a). Ultimately, they are combined into a single 
peak for L — > cx). Experimentally, the chain length of 
Sr2Cu03 is estimated ai L = 1800 ± 500, namely longer 
but finite. 
When the dimerization increases, the double peaks 



Fig. 20. System-size dependence of NMR intensity per spin (a) 
for the S = 1/2 Heisenberg model, and for CuGeOa (b) in the 
strong edge, and (c) in the weak edge. (S?)" is used. The impu- 
rity density is proportional to L~^ . The width of the gaussian 
is M) = 0.02. 



merge into a single peak, irrespective of the edge condi- 
tion. This single peak is attributed to the disappearance 
of delocalized staggered moment, because the intensity of 
the peak at i/ ^ originates from the dense distribution 
of {SfY near zero (see Fig. 13). There is one remarkable 
difference between the strong and weak edges; a broad 
skirt appears around the single peak for the latter, as 
seen in Fig. 19(b). This is caused by the localized free 
spins with large amplitude. 

Finally, we refer to the influence of the impurity con- 
centration. For (3=1 the similar aspects to the Heisen- 
berg case [see Fig. 20(a)] are observed for arbitrary values 
of a. In Fig. 20(b) we show the system-size dependence 
of g{H) for the spin-Peierls phase of CuGeOs. With in- 
creasing X {1/L) the height of the single peak lowers and 
the spectrum becomes broad. Such impurity effect seems 
common to various systems. 

§6. Haldane phase 

In this section, we study the edge effect in the Hal- 
dane phase, and clarify the difference from the spin liq- 
uid phase. As a typical compound, we refer to a F-AF 
bond- alternating compound monoclinic CuNb206. As 
discussed in Appendix D, this material is well described 
by the model eq. (1.1) with a ~ and f3 — —2.36 (—2) 
for the pure (10%-doped) sample. To avoid confusion. 



we make sure of the terms for edge condition: Corre- 
sponding to the strong (weak) edge in the spin-Peierls 
regime, we use terms "AF (F) edge" for (3 < 0, and "J 
(/3 J) edge" for a general value of /?. 

First of all, we consider the stability of the edge. In 
Fig. 21 we plot the energy difference AE in the ground 
states for a wide range of (3 and a = 0. As mentioned for 
the spin-Peierls regime, the strong edge becomes stable, 
as soon as f3 decreases from 1. However, as f3 further 
decreases beyond 0, AE starts to increase, and at last 
E[(3J) becomes lower than E{J) at /? ~ —3. For (3 -^ 
cx), AE becomes |/3|J/4 — oo, because the number of 
F bonds is larger by one in the F-edge case. Thus, in 
the S — \ Haldane systems the F (weak) edge is always 
realized. Incidentally, the whole feature of AE including 
the reversal point —3 only weakly depends on the value 
of a. In the bond-alternating systems with \(3\ <^ cxd 
like monoclinic CuNb206, in which the lattice structure 
is rigid in contrast to the spin-Peierls compounds, two 
kinds of edge J and [3 J should be yielded at random by 
non-magnetic impurities. 

Next, we discuss the spin polarization {Sf). As is 



Fig. 21. Energy difference between the J and (3 J edges in the 
ground states. To show little size dependence, some systems of 
even L are plotted together. E(,}) [E(I3J)] indicates the total 
energy for the system of the J (fij) edges. The minimum is 
AE/ J = —0.75 at /3 = 0. The arrow indicates the reversal point 
at 13 ~ —3. This value 3 is roughly the ratio of energy gain in 
singlet (3J/4) to triplet {J/4). 



known for the S' = 1 Heisenberg model, '^^'^*'' the mag- 
netic excitation is localized at both edges as S* = 1/2 free 
spins. This localized moment decays exponentially with 
the localization length ^odgc, which is independent of 
the system size and is equivalent to the bulk correlation 
length ^H ^ 6.03.^°'^®^ Magnetic LRO's hardly appear 
even by doping, unless the interchain couplings are suffi- 
ciently strong.^' ^'^' Besides, localized edge spins have 
been actually observed by various experimental stud- 
ies. '^^•' 

We have already seen that such localization appears 
also in the present model, when the dimerization is in- 
troduced for the weak edge (Fig. 13). Then, we start 
with the F (weak) edge. In Fig. 22 we compare {Sf) 



16 



Shinji Watanabe and Hisatoshi YOKOYAMA 



Fig. 19. Intensity of NMR spectra g{H) for a = (I) and 0.35 (II) for the segment of L = 200. The strong- and weak-edge cases 
are shown in (a) and (b), respectively. In 11(a) the line shapes are almost the same for /3 = 0.40-0.95. They are obtained from the 
staggered part of the data in Fig. 13, since (5^)" only weakly depends on i for /3 ~ 1. The width of the gaussian is ui = 0.02. 



in the Haldane phase (CuNb206) with the spin hquid 
(CuGeOs). For CuNbaOg (open circle), \{Sf)\ obeys 
typical exponential decay, in contrast to the spin-Peierls 
cases [see Fig. 13(b)]. Besides, ^edge is fairly short and is 
almost the same with the bulk correlation length ^ ^ 3.3 
(Fig. 5). By comparing Figs. 22(a) and (b), we find that 
the amplitude of (Sf) near the edge is independent of 
the system size. On the other hand, not only the oscil- 
lating part but also the uniform part of (Sf) vanishes 
for the bulk. Since the situation (Sf) = is the same 
with the ground state (5 — 0), the magnetic excitation 
never affects the bulk. Thus, the behaviors of the Hal- 
dane phase in the S = 1 systems stays unchanged for 
the present case. Incidentally, the odd-site systems show 
similar behaviors to this case, as discussed in §5.2. 

As for the AF (strong) edge, the excess spins extend 
to the whole system (small open diamond in Fig. 22). 
Nevertheless, the oscillating part is completely lacking. 
Thus, the magnetic excitations are not likely to induce 
LRO also in this case. Indeed, monoclinic CuNb2 06 
keeps the gap properties in heavily doped samples, ^^ in 
spite of the relatively large interchain coupling (cf. Ap- 
pendix D). In contrast, for CuGeOs the excited moment 
extends to the whole system (solid circle in Fig. 22). 
Moreover, the amplitude of staggered part abruptly 
grows with decreasing L (increasing x). Thus, the mag- 
netic excitation leads to AFLRO, as has been already 
discussed. 

Let us turn to the NMR spectra. In Fig. 23 we show 
the size dependence of the NMR intensity per spin for 
the parameters of pure CuNb2 06. For a low impurity 
concentration (large L) g{H) has a sharp peak at iJ ^ 
for every edge condition, corresponding to the absence 
of large oscillating part of (S'f ) in the bulk. For the F 
edge there appear skirts and broad side peaks due to the 
localized spins at the edges. Thus, the localization sur- 
vives even in the low-temperature limit, so that the NMR 



Fig. 22. Comparison of spin polarization in the lowest magnetic 
states between CuGeOs (strong edge) and CuNb2 06 (two kinds 
of edges) for the systems of (a) L = 100 and (b) 50. As for pa- 
rameter values, (o,/3)=(0.35,0.98) is assumed for CuGcOs, and 
(0.0,-2.36) for CuNb206. 



line shape will hardly broaden with decreasing temper- 
ature, unlike the S — 1/2 Heisenberg compound^*^ and 
the two-leg ladder system.^) On the other hand, when 
the impurity concentration increases (L becomes small) , 
the height of the single peak at H ^ lowers regardless 
of the edge condition. Simultaneously, the peak shifts 
slightly for the AF edge, while the intensity of the skirts 
and side peaks grow for the F edge. Anyway, the Haldane 
phase and the spin liquid share the tendency toward the 
breakdown of the sharp peak by a few percent of impu- 
rity. 

In conclusion, we recapitulate the relationship between 
the phase transition and the impurity effect. The char- 
acteristics of magnetic excitation in the Haldane phase 
are that the excited moment is confined in a finite range 
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Fig. 23. System-size dependence of the NMR intensity per spin 
for CuNb2 06 (a) in the AF edges and (b) in the F edges. For 
odd-site systems, g{H) is analogous to (b) with the doubled sys- 
tem size. We use the fixed value of /3 = —2.36 (for pure ma- 
terial) for simplicity, although the value of 13 changes with the 
impurity concentration (see Appendix D). To obtain the spectra 
by eq. (5.5), total (S^) (including uniform part) is used with 
w = 0.02. 



near the /3J edges, and that the oscillating part in (5f ) 
is almost absent for the J edges. The Haldane phase 
occupies the whole a-/3 plane except for the a axis and 
a = cx), on which the phase is the spin liquid charac- 
terized by the delocalization of (S'f ). As the parameters 
approach the spin liquid, the range of localized edge spins 
extends toward the bulk (independent of L) [(3 J edge], 
and the amplitude of oscillation (< 1/L) in (5f ) grows 
[J edge]. In the thermodynamic limit (for the pure ma- 
terial), however, both the range of edge spins [f3J edge] 
and the amplitude of staggered moment [J edge] are still 
negligible. On the other hand, when the chains are cut 
into finite segments (for the doped case), the range of 
edge spins may extend to the greater part of (or whole) 
segment [/3J edge] , and the amplitude of staggered mo- 
ment may becomes sufficiently large [J edge], according 
to the balance between (3 and L (x). Thus, the inter- 
chain coupling, which is quantitatively important, has a 
chance to induce AFLRO. In the transition areas from 
the Haldane phase to the spin liquid, impurities play a 
dramatic role. 

So far, we have not mentioned the two-leg ladder 
systems. In the typical ladder compounds, AFLRO is 
observed for the pure (LaCu02.5,*^°^ Tn ^ llOK) and 
hghtly-doped (SrCu203^'^^) cases. At the early stage 
these compounds were considered as isotropic ladders 



V^lcg/ ^r\ 



1), while recent experiments have 



concluded a ^ 2.^'^^' Actually, as described in §3, the 
isotropic ladder point [(a,/3) — (1,0)] exhibits a typi- 
cal properties of the Haldane phase with ^ '^ 3 without 
the tendency toward AFLRO. However, our preliminary 
calculations of various quantities as well as a previous 
study^**) indicate that the aspects for a ^ 2 arc rather 
similar to the S = 1/2 Heisenberg model. Thus, it is 
probable that the impurity effect of SrCu2 03 will be sim- 



ilarly understood by considering that the system is situ- 
ated in the transition area from the Haldane phase to the 
spin liquid in the two-chain Heisenberg limit {a -^ oo). 
We will publish the results for ladder elsewhere. 

§7. Summary 

We have studied a variety of properties of the ID 
Heisenberg model with the bond alternation and the 
NNN exchange coupling, using the DMRG and ED meth- 
ods. We summarize our main results below. As for the 
bulk properties before doping: 

(1) The conditions for the Haldane phase are satisfied 
on the Shastry-Suthcrland line (f3 — 2a) for < a < 0.5. 
A transition from the Haldane phase to a kind of spin 
liquid occurs at the Majumdar- Ghosh point (a — 0.5). 

(2) The above result can be extended to the whole 
a-(3 space. The Haldane phase continuously connects 
the S — I limit (/3 = — oo) to the spin-Peierls regime 
(/? '^ 1) through the bond-alternating regime, and the 
transition occurs at /3 = 1. In the Haldane phase, there 
is no anomaly in the dimer and string order parameters, 
whereas the spin gap and the correlation length exhibit 
singularities due to the incommensurate correlation. 

(3) We have determined the phase diagram of the in- 
commensurate correlation. The range of essential incom- 
mensurate correlation defined by (S^S^) is quite different 
from that by S{q). The singular points of the spin gap 
are included in the incommensurate range of Sij . 

(4) We have discussed various properties, especially 
dynamics, of monoclinic CuNb2 06. This compound is 
well described by eq. (1.1) as a typical ferromagnetic- 
antiferromagnetic bond-alternating chain. 

As for the edge effects, corresponding to the doping of 
non-magnetic impurities: 

(5) There are two kinds of chain ends, according as 
the edge bond is strong (J) or weak {(3 J). The strong 
edge is more stable for the spin-Peierls regime, while the 
weak (F) edge becomes stable for /3 ^ — 3. This fact is 
possibly connected to the compositional phase transition 
of the first order between the dimerized and uniform AF 
phases in CuGeOs. 

(6) The spin polarization is delocalized in the spin- 
Peierls systems and the two-leg ladders with large 
J\cg/ J-cnng, while the spin polarization is localized at the 
edge for the bond-alternating systems as well as the Hal- 
dane systems and the isolated dimers. In the former, 
(5*1) is enhanced by the impurity and triggers a 3D mag- 
netic long-range order. In the latter, {Sf) is insensitive 
to the impurity, so that the gap properties are preserved. 

(7) For finite but large systems NMR spectra in the 
low temperatures show a double-peak structure in the 
vicinity oi j3 ~ 1 due to the delocalized spin polarization. 
On the other hand, a single sharp peak appears for the 
Haldane phase, regardless of the edge condition. In this 
phase, the conspicuous broadening of the spectra with 
decreasing temperature will not be observed. Generally, 
as the impurity concentration is increased, the height of 
peak lowers and eventually vanishes. 
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Appendix A: Relations between two represen- 
tations 

In this Appendix, we summarize the relations between 
the Hamiltonians eq. (1.1) and eq. (2.3). Comparing 
corresponding terms, one obtains the relations. 



5 = 



1-/3 



2a 



J = 



1 + /3 



J. 



(Al) 



1 + 13' 1 + /3' 

Thereby, a line (point) in the a-/? plane is moved to 
another line (point) in the a-d plane. Note that the 
Shastry-Sutherland line is /? = 2a in eq. (1.1) and 
2a + (5 = 1 in eq. (2.3). The isotropic ladder point is 
(a,/3) = (1,0) or {a,S) = (2,1). However, the a axis 
(a,/3 = 1) remains the a axis {a,S = 0). The energy 
scale also should be changed; E in eq. (1.1) is related to 
E in eq. (2.3) as. 



E 
J 



2 E 

1 + 137' 



(A-2) 



Appendix B: Cusp Behavior of the Gap 

In this Appendix we pursue the cusp behavior in A to 
supplement the discussions in §3.1. 

The second-order perturbation energies for the ground 
state (5 — 0) and the lowest excited branches {S = 1) 
around the isolated dimcr point are given by eqs. (3.1) 
and (3.2), respectively. The singlet energy is lowered by 
the second-order perturbation for (3 ^ 2a; but this decre- 
ment vanishes for f3 — 2a (the SS line). This cancella- 
tion is understood intuitively as follows: Two a bonds 
and one /? bond connect two singlet dimers, but a and /3 
play mutually inverse roles for the spin directions in the 
two singlet pairs. On the other hand, the triplet branch 
has a dispersion due to the motion in the array of singlets 
for P 7^ 2a. The bottom of Et{q) is located a.t q ~ for 
P > 2a and a,t q — n/2 for /3 < 2a; the IC value of q 
does not appear around the isolated dimer point. Thus, 



A is obtained as: 
1 



A 

7 



1 



1 



f 

-(/3-2a) 

1 ^^ 

-(/3-2a) 
3 



3 2 3 

-a — -pa 

4 4 



/3 = 2a 



(3>2a 



lb A A 

Comparing the first-order terms, one finds the maximum 
of A is at /3 = 2a. The cusp behaviors including the 
curvature in both sides oi (3 — 2a are reproduced by the 
differential coefficients: 




a^A 



<0 



>0 



I3>2a 
(3 < 2a. 



With increasing a, however, the cusp point deviates to 
the side of /? < 2a as seen in Fig 12. 

Then, let us look at the cases of finite a with the ED 
results. In Fig. 24(a) the lowest triplet excitation for 
the available momenta in L — 24: are plotted for a — 
0.35. The boundary between the lowest excitations of 
q = TT and q = tt/2 still forms the cusp near /3 = 0.525. 
An interesting point is that the excitation energies for 
all the momenta meet at the cusp just like the isolated 
dimers. No dispersion in the lowest excitation is evident 
in S{q,u!) (Fig. 25); the excitation has a local character 
at the cusp points. 

The above feature continues to a ~ 0.5, beyond which 
the lowest excitation is no longer restricted to the C wave 
numbers. As in Fig. 24(b) for a — 0.6, the wave number 
of the lowest excitation docs not switch directly from 
q = 7r/2 to TT, but IC excitations lie between them. Thus, 
the cusp behavior disappears for a ^ 0.5, as seen for 
a = 1 in Fig. 3. 



Fig. 24. Excitation energy of the lowest triplet branches for each 
available value of q for L = 24; (a) a = 0.35 and (b) 0.6. The 
data are obtained by ED and the recursion method. 
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Fig. 25. Dynamical structure factor S{q,iLi) for a cusp point 
{a, /9)=(0. 35, 0.525). Intensity is proportional to the area of cir- 
cle. Data for L = 14-26 are plotted together. Note that every 
branch in this figure is independent of L. 



Appendix C: Low-Energy States on the Shastry- 
Sutherland Line 

The isolated dimcr point satisfies the conditions of the 
Haldane phase itemized in §2, while the MG point breaks 
the condition of a unique ground state. Hence, the phase 
is different between these points, which are the two ends 
of the SS line. With a phase transition in mind, we 
summarize and reconsider the knowledge of the ground 
state and low-lying excitations along the SS line in this 
Appendix. In §C.l using a non-local unitary transfor- 
mation, we explain that the SS line strictly satisfies the 
conditions of the Haldane phase. In §C.2 we summarize 
the structure of energy levels under different boundary 
conditions with the ED data. In §C.3 we show that the 
transition occurs just at the MG point in terms of the 
perturbation theory for the degenerate case. 

C.l Non-local unitary transformation 

Following the formulation of the non-local unitary 
transformation by Takada and Kubo,^^'^^^ we rewrite 
the Hamiltonian eq. (1.1) as 

L-l 



= 1 
L-l 

E\ .^X —X I .— Z __Z I __3^_.iC I -—Z^_,Z 



i=l 
L-l _ 

a 






aY.^. 



(C-l) 



where the operator S is replaced by CTi = S2i and Ti — 
S2i-i- As for the isolated dimer point (a = /3 = 0), the 
surviving part [the first line of eq. (C-l)] is coupled Ising 
chains with a Z2 x Z2 symmetry. The ground state has 
ferromagnetic orders in both erf and rf , and is fourfold 
degenerate {af = ±1/2 and rf = ±1/2); the symmetry 



is completely broken.^*' ^^^ Moreover, since the eigenstate 
is V'l of eq. (2.1) (or 1/12), there exists a single-triplet gap 
of magnitude J. The string order parameter is calculated 
as 1/4. Thus, the conditions of the Haldane phase are 
satisfied at this point. These conditions are still satisfied 
on the SS line for < a < 0.5, because the ground state 
is identical with that of the isolated dimer point, and is 
unique with a finite singlet-triplet gap. We will make 
sure of the latter point below. This result has been also 
reached with the matrix product procedure. ^^ 

C.2 Structure of low- energy levels 

Here, we study the low-energy levels on the SS line 
((5 + 2(5 = 1) for finite systems. The boundary condition 
greatly affects the level structure. 

Figure 26(a) shows the ED result under PBC in the 
S^ = sector. Under this boundary condition there is 
no degeneracy in the ground state in accordance with the 
conditions of the Haldane phase. At the isolated dimer 
point, the first excited states are triplet and 3L/2-fold 
degenerate. Although this degeneracy is lifted by a fi- 
nite value of a, the energy difference is very small except 
for the vicinity of the MG point. There, as seen in the 
inset, the level of S* = 0, which is one of the degenerate 
ground states at the MG point, goes up linearly with 
decreasing a, and intersects the lowest S — I level at 
a ~ 0.48 (indicated by arrow). This means that the 
lowest excited state becomes singlet near the MG point. 
However, the range of this singlet-singlet excitation de- 
pends on the system size. In Fig. 27 we show similar 
figures with different system sizes. With increasing L, 
the crossing point approaches the MG point. This point 
will be pursued in §C.3. 

In Fig. 26(b) plotted is the result for OBC with the 
strong edges. In this case the ground state is either ^1 
or i/)2 according as singlet pairs sit on the edge bonds, and 
always nondegenerate. Note that the edge bonds are oc- 
cupied by singlet pairs even for the non-dimerized (MG) 
case. The first excited state is always triplet. Thus, not 
only the ground state is identical with that of PBC, but 
properties of excitation are also analogous to those of 
PBC. We have taken advantage of such a property to 
obtain the bulk gap in §3.1. 

Figure 26(c) shows the same for OBC with the weak 
edges. In this case, there is the pseudo fourfold (twofold 
in the S'^ — sector) degeneracy except for the vicinity 
of the MG point. The ground state {S = 0) is neither tpi 
nor V'2 except for the MG point, due to the isolated spins 
at the edges. Thus, this edge condition has a feature spe- 
cific to the open edge of the Haldane systems. The nearly 
degenerate lowest levels (S* = and S — I) separate at 
a ^ 0.46; the S ~ level abruptly comes down and con- 
verges to the value of V'l (or V'2) at the MG point. In this 
process the singlet array is rearranged, say from L to R, 
to adapt itself to the strong edges. On the other hand, 
near the above separation point the second lowest S = 
level, which has a similar energy to the lowest S — 1 
level near the MG point, abruptly changes its direction 
and merges into the upper S = I levels. These behaviors 
correspond to the level crossing under PBC mentioned 
above, and the separation point also highly depends on 
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Fig. 26. Some lowest levels in the S^ = subspace on the SS line 
(<5 + 2o = 1) as a function of a; (a) PBC (L = 16), (b) OBC with 
the strong edges (L = 12), and (c) OBC with the weak edges 
{L = 12). The insets are respective magnifications near the MG 
point. The arrow in the inset of (a) indicates the crossing point 
of the second lowest 5 = level and the lowest 5 = 1 level. 



Fig. 27. Three lowest levels on the SS line (5 + 2a = 1) near the 
MG point in the 5"^ = sector under PBC. The system size is 
different between (a) and (b) [see also the inset of Fig. 26 (a)]. 
The Arrow indicates the crossing point of the lowest 5 = 1 level 
and the second lowest 5 = level. 



the system size. This point will be taken up again in 
§C.3. 

In this connection, recent studies*^^ have discussed the 
appearance of the mid-gap states for the doped CuGeOa 
at Tn < T < Tgp. According to the above results, how- 
ever, mid-gap levels are not likely to arise from the strong 
edges, the case for CuGeOs. On the other hand, for the 
weak edges (including odd-site systems) the energy dif- 
ference among the nearly fourfold degenerate states in- 
creases exponentially with l/L. Thus, the experiment 
may detect such levels, if any. 

Finally, we generalize the level structure for the bond- 
alternating systems, aside from the behaviors near the 
MG point. Under PBC or OBC with the strong edges, 
the ground state is a nondegenerate singlet (^i or ip2 
in the case of the SS line) with the singlet-triplet gap 
over it. Under OBC with the weak edges, the ground 
state is nearly fourfold degenerate (one singlet and three 
triplets) with the bulk gap over these four states. 

C.S Breakdown of the perturbation theory 

Having seen extraordinary behaviors around the MG 
point in §C.2, here we reconsider them in the light of 
the perturbation theory for the degenerate case, and 
shows that the non-Haldane phase is restricted to the 
MG point. 

The Hamiltonian eq. (2.3) on the SS line is written 
around the MG point as. 



7i = Hmg + CH', 



L r 



Wmg = J 7 ^ 



i=l 



Si ■ Si^i -^ —Si ■ Si^2 



L r 



w 



-2^E 



i=l 



(— l)*^^ • Si+i + —Si ■ Si+i 



(C-2) 



(C-3) 



, (C-4) 



where we assume PBC and C = 1/2 — a to be infinites- 
imal. Here, TLyiG'^n = ^McV'n {n = 1,2) holds, where 
Emg — —3LJ/8 is the ground-state energy at the MG 
point, and ipn is assumed to be normalized. Let ipi be 
the ground state on the SS line: Ti'ipi = eiV'i with 
El = — 3ivJ/4.^^^ Since -01 and ■!/'2 are orthogonal in the 
thermodynamic limit and ipi is still an eigenstate under 
the perturbed Hamiltonian, the first-order perturbation 
energies for the two levels are formally written in a sim- 
ple form, as 



K = C(V'n|H'l^n)=Cen (^-1,2). 



(C-5) 



Thus, the perturbation theory naturally works for the 
lower level ip[ — ipi- As for the higher level ■02) if the 
perturbation theory works, the energy is also propor- 
tional to 1/2 — a near the MG point, actually as in the 
inset of Fig. 26(a) and in Fig. 27. Since a direct calcula- 
tion of £2 looks difficult, we estimate it by ED with PBC 
through. 



£2 



(V'2|W>2) 



dE. 



dC 



(C-6) 



C--0 



where E2 is the total energy with respect to tp2 . 
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In Fig. 28(a) we plot J/e2 thus obtained versus \/L. 
It shows that £2 diverges in the thermodynamic hmit, 
namely the perturbation breaks down for ^2 • Along with 
this phenomenon, the crossing point of ^2 ^^id the low- 
est S = \ level converges to the MG point [Fig. 28(b)]. 
Thus, one of the degenerate levels at the MG point never 
continues to the dimerized region [5 > 0); the first ex- 
cited state under PBC is always triplet on the SS line. 

On the other hand, as pointed out for OBG with weak 
edge in §G.2, the lowest level, which separates from the 
triplet levels near the MG point, behaves similarly with 
V'2 in PBC. We also estimate the perturbation energy e 
for this case in the same manner, and plot J /e together 
in Fig. 28(a). Within the numerical accuracy, the val- 
ues of J/e2 and J/e are identical. Thus, the lowest level 
(5* = 0) at the MG point does not smoothly continue to 
the 5* = state which makes up the fourfold degeneracy 
for (5 > 0. And the degeneracy of the ground state con- 
tinues to the limit a -^ 1/2. Thus, the energy difference 
between the strong and the weak edges in the limit of 
L ^ 00 and a — > 1/2 becomes the MG gap. 

From the above we conclude that the Haldane phase 



Fig. 28. (a) System-size dependence of the inverse of the first- 
order perturbation energy per site for ^'2 i^i PBC. The corre- 
sponding value for the lowest level in OBC with the weak edges 
are also shown, (b) System-size dependence of the crossing point 
of 1P2 and the lowest S = 1 level in PBC. 



continues to this limit, and that the MG point is critical. 
Hence, the consequence of §C.l is valid for S ^ 0. On the 
other hand, we can say inversely that in finite systems 
the properties of the Haldane phase are lost near the MG 
point. So far, we have restricted our discussion to the SS 
line. As far as we have tried numerically, however, this 
discussion is applicable to any path which starts from an 
arbitrary point on the a axis (a > etc) and advances in 
an arbitrary direction of finite S. 

Appendix D: Bond-alternating compound 
CuNb206 

In this Appendix we discuss various aspects of a F-AF 
bond-alternating chain monoclinic CuNb2 0e to supple- 
ment the background of §6. The dynamical properties of 
bond-alternating chain, in particular for /3 < 0, are sur- 
veyed in §D.l. In §D.2 various properties of CuNb206 
are summarized, before discussing its dynamical proper- 



ties in §D.3. 

D. 1 Spin dynamics for a ~ 

Since principal static properties have been already de- 
scribed in this paper as well as in some studies,^^'^®'*'*) 
we pay attention to the dynamics. As discussed in §5.1, 
low-lying spin excitation change its character from even- 
number spinous to magnons, as soon as the bond al- 
ternation is introduced. Consequently, isolated magnon 
branches appear below the continuum. As in Figs. 29(a) 
and (b), the ratio of intensity (residue) in the lowest iso- 
lated branch increases with decreasing (3 from 1 . At the 
isolated dimer point this branch has the total intensity 
with a completely fiat band at a; = J. For /3 > the 
wave number of gap is tt and the band top is at tt/2. 

On the other hand, as /3 decreases further for /3 < 0, the 
main branch comes to have dispersion with a minimum 
at q ~ 7r/2, as discussed by the perturbation calculations 
in Appendix B. A continuum due to two triplet pairs 
appears around lu = 2J, as pointed out for (3 ^ by 
Hida.^^-* Furthermore, another isolated branch appears 
approximately at a; = |/3|J, which corresponds to the 
energy to break one F bond. In addition, higher har- 
monics appear, e.g. at uj ^ 6.5 J in Fig. 29 (d). 

To begin with, let us look at the dominant lowest- 
energy branch (LEB) and the continuum. As shown in 
Fig. 10, the maximum of integrated intensity S{q) moves 
from q = TT to -k 12 through IC wave numbers, as /3(< 0) 
decreases. According to this change, and because the 
greater part of intensity is concentrated on LEB, the 
weight in LEB moves from (7 = tt to 7r/2. For /3 ^ — 3 the 
intensity at g ~ 7r/2 is the strongest, and grows further 
with decreasing /?(< 0) [Fig. 30(a)]. Similarly, the weight 
of the continuum chiefly distributed at g ^ tt for small 
values of |/3| gradually moves to g ^ 7r/2. We have found 
for /3 ^ - 2 that the residue of LEB for 0.2 < qj-K < 0.5 
is more than 85% and scarcely depend on the system 
size, whereas for q ^ 0.27r it considerably decreases with 
increasing L. This implies that LEB enters the contin- 
uum at q ^ 0.27r. A similar aspect was observed for the 
S = 1 Heisenberg model. ^^^ 

Next, we mention the high-energy branch (HEB) near 
ui ^ \l3\J. The dispersion of HEB is flatter with a mini- 
mum also at q/2. Accordingly, local destruction of a fer- 
romagnetic bond chiefly contributes to this excitation. 
There exist weak continuum spectra around HEB, al- 
though HEB itself is an isolated branch because of little 
dependence of the intensity on L. As shown in Fig. 30(a) 
the intensity of HEB increases with decreasing /3 at first; 
the residue for HEB at q = n becomes as large as 70% 
for /3 < —4 in spite of the high-energy process. With de- 
creasing /3(< 0) further, however, its intensity decreases 
and vanishes in the Haldane limit, according to the to- 
tal intensity 5(7r) (Fig. 10). Therefore, this branch is 
specific to the F-AF bond-alternating systems. 

D.2 Monoclinic CuNb2 0e 

As a spin-gap compound CuNb206 has been ac- 
tively studied by various experimental means re- 
cently.^' ^^'*^"*^-' CuNb206 has zigzag chains of edge- 
sharing CuOg octahedra with localized S* = 1/2 spins 
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Fig. 29. Dynamical structure factor S{q, ui) for four values of /B (5) and a = 0. Data obtained by the ED with a recursion method 
for L = 14-26 are simultaneously plotted. The intensity is proportional to the area of circle. In (a) and (b), to convert the a-S 
representation into the o-/9 representation, 5/6 and 2/3 have to be multiplied to w, respectively. See Appendix A. 



Fig. 30. (a) Intensity [S{q, lu)] of LEB at q = 7r/2 (solid circle) 
and of HEB at 5 = tt (open circle) as a function of — /3. The data 
points are extrapolated from the ED results of L = 14-26. (b) 
Excitation energy of LEB at q = 7r/2 [u}i,{= A)] and of HEB at 
q = n [t<JH]) and their ratio are shown as a function of — /3. These 
values are estimated by the DMRG and ED methods (L -^ 00) 
for eq. (1.1) with o = 0. The experimental value of the ratio for 
10%-doped CuNb2 06 are indicated by the dashed line. 



on Cu^+ ions. This compound has two polymorphs, i.e. 
orthorhombic and monochnic phases. The orthorhombic 
phase is considered as a ID Hcisenbcrg antiferromag- 
net with an AFLRO below 7.3K; on the other hand, the 
monoclinic phase shows thermal-activation-type behav- 
iors in magnetic susceptibility x{T)-^^' Henceforth, we 
restrict our discussions to the monoclinic phase. Ow- 
ing to the zigzag structure and the neutron scatter- 
ing data of the lowest excitation at q — tt/2, it is 
expected that monoclinic CuNb2 06 should be a F-AF 
bond-alternating chain. Actually, the experimental xC^) 



was fitted successfully to the results of the model eq. (1.1) 
with J = 51K, a = and (3 = —2.36.*^ Simultaneously, 
the coincidence in specific heat C{T) is satisfactory with 
the same parameter set. Our estimation A/ J — 0.3714 
for these values (Fig. 3) leads to A = 19K. 

Notwithstanding, recalling the ambiguity in modeling 
a spin-gap systems based solely upon x(2^)j for example 
(VO)2P207*^) and KCuClg,^") one should run a careful 
check by other means, in particular neutron scattering. 

Before discussing the neutron experiment, we touch 
on a couple of related properties of CuNb206. First, 
we consider the effect of the NNN coupling. Fitting 
the experimental xi^) to the ED results of L = 16^^^ 
around (a,/3) = (0, —2.36), we have found that a series 
of parameter sets such as (0.10,-2.0), (0.24,-1.5) and 
(0.38, —1.0) bring favorable results in the same degree 
with (0,-2.36). This reminds us of the tendency that 
a encourages the TTii correlation. Concerning C{T), 
however, such a set leads to a poorer fit, as a increases. 
Hence, we conclude that the NNN exchange is negligibly 
small in CuNb2 06 despite the edge-sharing linkage. This 
is intuitively understood by considering the difference in 
the orbital structure between CuGeOa and CuNb2 06- In 
CuGeOs the NNN exchange arises from a superexchange 
process through the path of Cu-0-O-Cu in the Cu02 
plane. ^^^ In this case a hybridized Cu-0 orbital overlaps 
the counter Cu-0 orbital to some extent. On the other 
hand, in CuNb206 the two Cu-0 hybridized orbitals of 
our concern are not in a plane and do not overlap each 
other because of the zigzag linkage of CuOg octahedra. 
This comparison affords vital support to the above mech- 
anism of the NNN exchange in the edge-sharing com- 
pounds. 

Next, we consider the effect of elemental substitu- 
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tion.*' In the crystal of monoclinic CuNb206, Cu can be 
replaced by Zn up to 40%. In sharp contrast to CuGeOs, 
the gap behaviors remain and no magnetic LRO ap- 
pears even for the 40%-doped sample. In doped sam- 
ples, however, lattice constants and bond angles change 
according to the Zn concentration, mainly due to the 
different atomic radii between Cu and Zn; this change, 
particularly in the bond angles, ought to affect the ex- 
change integral, according to the Kanamori-Goodenough 
rule.^^-* Actually, considerable changes in parameters oc- 
cur as summarized in Table 11,^'^) where the values are 
estimated by the fit of x{T) for some doped samples 
Cui_a;Zn2;Nb206. Also entered in Table II is the gap es- 
timated through eq. (1.1) using those parameters. The 
decreasing tendency versus x coincides with the gap es- 
timated from C{T), which decreases linearly with x.^' 
This behavior is contrary to what the S — 1 Haldane sys- 
tems show; the gap increases with increasing x.^^' The 
reason is that thick doping cuts the chains into shorter 
segments, and that the gap generally increases as the 
system size becomes small. In the case of CuNb2 06 the 
decrease of | J|, which is specific to this compound, prob- 
ably overcomes the system-size dependence of the gap 
and the increase of /3. 



Table II. Parameters for the doped CuNb2 06 obtained from the 
fit of xC^)- The gap A is estimated through the model eq. (1.1). 



X 


J(K) 


/3 


A(K) 


0.0 


51.4 


-2.36 


19.1 


0.1 


44.5 


-1.95 


18.9 


0.2 


40.0 


-1.86 


17.6 


0.3 


33.4 


-1.43 


17.2 


0.4 


28.1 


-1.03 


17.1 



D.S Inelastic neutron scattering 

Recently an inelastic neutron scattering experiment 
was performed on an Zn-doped CuNb2 06 sample. ^^•' The 
spectra obtained in this experiment have reproduced the 
most features of the F-AF bond-alternating chain dis- 
cussed in §D.l. We owe all the experimental data in this 
subsection to ref. 15. 

We begin with the evaluation of the parameters. In 
Fig. 30(b) shown are the excitation energies of LEB wl 
at q — tt/2 (o^l — A), and of HEB o^h at q = tt, and their 
ratio wh/wl, estimated by the DMRG and ED methods 
for eq. (1.1) with a = 0. On the other hand, the experi- 
mental values are lu' 



L 



1.73meV and uj^ 



13.14meV, 

respectively. Then, the ratio becomes uj^/uif^ = 7.6, 
which is shown with dashed line in Fig. 30(b). Thus, we 
find that the value /3 for CuNb206 of 10% substitution is 
about —2.0, which agrees well with the values estimated 
by x{T) in Table II. 

In Figs. 31(a) and (b), we depict the excitation energy 
of LEB and of HEB, respectively, for the whole q space. 
In these figures wc also plot the experimental results with 
J = 46. 9K, which best fits the whole data (both for 
LEB and HEB) at the same time to the calculated values 



for (3 = —2.0. They agree well as a whole except for 
q — O.Stt in LEB, where the experimental accuracy is 
rather poor due to the weak intensity. The above value 
of J is consistent with J = 47. 9K obtained from the 
gap A. In Fig. 31(b) we also show the data of a ^ 
and /3 = —2.0. Although HEB is concave for a = 
irrespective of the value of /3, as seen in Figs. 29(c) and 
(d), yet it becomes convex as a increases. This feature 
for finite a contradicts the experimental results. Thus, 
the neutron experiment supports the absence of the NNN 
exchange. 
In Figs. 32(a) and (b) shown are the intensity [5'(g,a>)] 



Fig. 31. Excitation energies (a) of LEB for some values of /3, and 
(b) of HEB for a few values of /3 (o = 0) and for some values of 
a (/3 = —2.0), as a function of q. The data for L = 14-26 are 
simultaneously plotted, which are obtained by ED and recursion 
methods. A large gray circle with a error bar represents the 
experimental data of 10%-doped CuNb206^'^' with J = 46.9K. 



of LEB and of HEB, respectively. In these figures we 
scale the experimental results so that LEB may be best 
fitted. For LEB in Fig. 32(a) the experimental inten- 
sity is roughly reproduced by the calculated values for 
P = —2.0. Kodama et al. pointed out the conspicuous 
intensity at q — 0.67r.^^^ This is probably because the 
system (/3 = —2.0) is in the regime of IC correlation 
in S{q) as discussed in §4.2. Then, S{q) has a peak at 
q = O.SIStt (See Fig. 11). On the other hand, the exper- 
imental intensity for HEB [Figs. 32(b)] is considerably 
smaller than what the theory expects. This discordance 
is still an open question. 
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Fig. 32. Intensity [S{q, lu)] of (a) LEB for several values of /9, and 
of (b) HEB for a few values of /3; both for a = and as a func- 
tion of q. The data for L = 14-26 are simultaneously plotted, 
obtained from the ED and recursion methods. The diffuse data 
points for /3 = —1.5 in (b) is owing to the overlap with the con- 
tinuum spectrum around HEB. A large gray circle with a error 
bar represents the experimental value of 10%-doped CuNb2 06 
scaled so as to best fit the lowest-energy branch of /3 = —2.0. 
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